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ABSTRACT 


Probabilistic composite mi ciu mechanics methods are developed that 
simulate expected uncertainties in unidirectional fiber composite 
properties. These methods are in the form of computational procedures 
using Ifcnte Carlo simulation. The variables in which uncertainties are 
accounted for incline constituent and void volume ratios, constituent 
elastic properties and strengths, and fiber misalignment. A 
graphite/epoxy unidirectional composite (ply) is stulied to demonstrate 
fiber composite material property variations indited by random changes 
expected at the material micro level. Regression results are presented 
to show the relative correlation between predictor and response 
variables in the study. These computational procedures make possible a 
formal description of anticipated random processes at the intraply 
level, arai the related effects of these on composite properties. 
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CHAPTER I 


IHTRODUCTION 


A. Background 

The diverse requirements of recent engineering applications have 
motivated designers to explore specialised structural and material 
systems. Ceramic naterials, for exanple, have several attractive 
structural properties, such as their high stif fnessAcight ratios, and 
low variation of stiffness and strength over wide ranges of 
environmental conditions. A significant disadvantage inherent to 
brittle structural naterials is their vulnerability to failure die to 
cracks propagating from flaws. The increased probability of a flaw 
occurring in a material as the volune increases leads to bulk strengths 
vfcich are a fraction of the tteoretical strength of the material. The 
size effect on naterial strength (Ref. 1) can be explained by the 
-weakest link" concept. Griffith (Ref. 2) reasoned that very small 
solids, for exanple wires or fibers, might be expected to be stronger 
than large ones, due to the additional restriction on the size of the 
flaws. In the limit, a single line of molecules mist possess the 
theoretical molecular tensile strength of a material. A conseqimnce of 
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the size effect on strength was the development of fiber composite 
materials tfiich consist of thin, strong fibers bound together by a 
d«=tile matrix. The advantages of fine, strong fibers can explain the 
current trend toward increased use of fiber composite materials in 
demanding aerospace applications. 

Properties of a composite laminate depend on the properties of the 
constituent materials, their distribution, and orientation. Laminates 
are composed of layers of unidirectional ly reinforced plies (laminae). 

The lamina is typically considered the basic unit of material in a 
composite strixjtural analysis, vfcich requires knowledge of the material 
properties of each individual lamina and its geometric orientation. The 
branch of composite mechanics that predicts ply material properties 
based on the properties, concentration, and orientation of its 
constituents is known as composite micromechanics , and frequently 
irexmporates the traditional ffechanics of Materials assumptions. The 
desired laminate is created by stacking of plies in specific directions. 
The integration of ply properties to yield laminate properties is called 
laminate theory. Laminate variables six* as ply orientation and 
stacking seqi*nce can be tailored to yield a laminate with the desired 
material properties. Thus, the laminated composite is a suitable 
material for component design- 

Analysis of fiber composite structures is currently performed using 
a variety of computer codes. From the original codes based on classical 
micromechanics and laminate theory, recent codes (Bef. 3,4) have been 
developed idiich incorporate the current state of the art. Complete 
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mechanical, thermal, and hygral properties are calculated, and can be 
used to compute response. Advanced failure criteria are used to 
calculate composite strengths. Environmental effects are also 
quantified. The usefulness of these codes has been demonstrated by 
comparison with experimental and finite element results (Ref. 5,6). 

Tte analytical capability of many codes is limited by the 
deterministic nature of the computations. Specifically, fixed values 
for constituent material properties, fabrication process variables (i.e. 
constituent volume ratios) and internal geometry mist be used as input. 
Ifeuever, random variations in these parameters are not -only expected, 
but easily observed experimentally. (See Fig. 1) 

The analysis of composite structures requires reliable predictive 
mcxlels for material properties and strengths. Hovever, the prediction 
efforts have been complicated by inherent scatter in experimental data. 
Since vreertainties in the constituent properties, fabrication 
variables, arxl internal geometry vould lead to uncertainties in the 
measured conposite properties, the qiestion arises: 

How much of the "statistical*' scatter of experimentally observed 
composite properties can be explained by reasonable statistical 
distribution of input parameters in composite micronechanics and 
laminate theory predictive models? 

The increasing use of probabilistic methods in structural mechanics has 
been shewn to provide a more realistic depiction of structural response 
due to load variations. (Ref. 7) The recognition that material 
parameters are characterised by a spectra of values ( that is, are 
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statistical in nature ) rather than by a unique set of values, points to 
probabilistic nethods as a logical analysis approach. 



s 



Fig. 1- 


Photonlcrograph of Graphite/Epoxy cross section 

showing variation in fiber content. IRef. i?) 
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B. Purpose 

The aim of this thesis is to develop a computational oapability to 
simulate the probabilistic variations in the mechanical behavior of 
wudirectional fiber composites. The Ifonte Carlo method is used to 
simxlate a variety of random processes , to qtantify fiber cooposite 
material variations induced by random changes in composite fiber 
alignment, constituent properties, and fabrication process variables. 
This ranlom process description is an attempt to more accurately predict 
thp behavior of manufactured materials, %hich inherently inclufe these 
ranlom variations. The characterization of fiber reinforced composites 
through simulation of local nonuni formi ties provides an economical 
alternative to experimentation to measure material properties. 
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C. Formulation of the Model 

The nodel oomumly used in characterizing fiber conposites is basmi 
„ the calculation cf properties of the basic unit of an orthotropic 
ply. The layup georetry is then used in laminate equations to calculate 
conposite properites (See Figs. 2a, 2b). In this -rh, ho«uer, the 
basic unit is taken as the sub-ply, shinh consists of only one 
fiber-retriH level in the material. flicronechanics theory is used to 
calculate the properties of the assured orthotropic sub-ply. each with 
randomly distributed fabrication variables and material properties. 
Distributed fiber directions, dre to possible misalignrent within the 
ply, are then used in the laminate equations to calculate ply 
properties, this substrreturing cf the conposite ply represents a novel 
attenpt at characterization of fiber conposite reterial properties based 
on probabilistically distributed constitrent properties, individual 
fib«- misalignrent, and fabrication process variables (See Figs. 3a, 3b). 

This formation is particularly well suited to the probabilistic 
description of fiber conposite reterial properties. Since the 
micronechanics and laminate equations ere be used to calculate ply 
properties at any uunher of points in a ply, a tractable finite elerent 
structural analysis based only on sinple distributional assunptiens tor 
physical parareter variations can be per f creed. This nodel supplies a 
rational procedure for conposite reterial property assessnent, 
it treats the reterial as the result of a series of random processes 
at the intraply level. 


which occur 
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(a) orthotropic ply (b) laminate 

Fig. 2- Conventional Model 



(a) subply 



Fig. 3- Substructure Model 
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D. lfethod of Investigation 

1 . fk'ief Description of ICAN 

The Integrated Composite Analyzer (ICAN) is a computer program for 
comprehensive linear analysis of multilevel fiber composite strictures. 
The program contains the essential features required to effectively 
design structural components made from fiber composites. It now 
represents the culmination of research conducted since the early 1970 s, 
at the National Aeronautics and Space Administration (NASA) Lewis 
Research Center (LeRC), to develop and code reliable composite mechanics 
theories. This user friendly, publicly available code incorporates 
theories for 


1. con v entional laminate analysis 

2. intraply and interply hybrid composites 

3. hygral, ttermal, mechanical properties and response 

4. ply stress-strain influence coefficients 

5. microstresses and microstress influence coefficien s 

6. stress concentration factors around a circular hole 

7. predictions of delamination locations around a circular hole 

8. Poisson’s ratio mismatch details near a straight free edge 
g. free edge interlaminar stresses 

10. laminate failure stresses 

11. normal and transverse shear stresses 

12 explicit specification of matrix-rich interply layers 
£ finite eleS^t material cards for NASTRAN, HARC 


A detailed description of ICAN can be found in Reference (3). The 
ICAN code and documentation are available through COSHIC, the Computer 
Software Ifenagement and Information Center, Suite 112, Barrow Hall, 

Athens GA , 30602 - 
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2. Sumnary of Variables 

The variables studied in this work can be separated into two 
categories. The independent variables to be simulated using random 


sampling consist of the following (see 

Fig. 4a for fiber coordinate 

system) : 


Geometry: 

(THETA) 

fiber orientation angle 

Fabrication variables: 


fiber volune ratio 

(FVR) 

void volune ratio 

(WR) 

Fiber properties 

(EFP1) 

longitudinal elastic modulus 

transverse elastic modulus 

(EFP2) 

shear modulus, 1-2 plane 

(GFP12) 

shear modulus, 2-3 plane 

(GFP23) 

fiber tensile strength 

(SFPT) 

fiber compressive strength 

(SFPC) 

Matrix properties 

(EM*) 

elastic modulus 

matrix tensile strength 

(SIFT) 

matrix compress ite strength 

(smpc) 

matrix shear strength 

(SIPS) 

Tte deperxlent variables to be calculated using I CAN consist of 

following ply properties, measured about the material axes (see Fig 

4b): 


normal modulus in 1-1 direction 

(ecu) 

normal modulus in 2-2 direction 

(EC22) 

shear modulus in 1-2 plane 

(EC12) 

Poisson's ratio for strains in 2 direction induced 

by stresses in 1 direction 

(N0C12) 

Poisson's ratio for strains in 1 direction induced 

by stresses in 2 direction 
Coefficients of thermal expansion 

(NUC21) 

in 1-1 direction 

(CTE11) 

in 2-2 direction 

(CTE22) 

coupling coefficient 

(CTE12) 
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SlCAJt 

(a) fiber (b) material 


Fig. 4- Coordinate Systems 



Fig. 5- Order of ICAN input data cards 
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Ply strengths in material directions 
longitudinal tensile 
longitudinal cowpressiw 
transverse tensile 
transverse ccnpressiw 
in-plane shear 


(9CXXT) 

(SCXXC) 

(scyyt) 

(SCYYC) 

(SCOTS) 


descriptions above should be consulted periodically for the 
definitions of variables that henceforth will be referred to 
syntolically. 
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3.fbnte Carlo Methods 

Complicated stochastic processes can be si mil a ted by a variety of 
nunerical methods generally referred to as ffante Carlo aethods (Ref. 8). 
Hie term refers to that branch of experinental mathematics concerned 
with experiments on random nunbers. Since the adrent of high speed 
computers, they have found extensive use in most fields of science and 
engineer ir^, in analysing many physical processes of a statistical 
nature, or vtere direct experimentation is not feasible. In general, 
they can be economically used to achieve a level of precision between 90 


and 95 percent . 

ft rfonte Carlo experinent refers to the procedure of randomly 
assigning a value to an independent random variable in a chosen model , 
and observing the dependent variable at the conclusion of the process 
being modeled, ft Monte Carlo procedure is composed of n stch 
independent experiments. Wien n is sufficiently large, the observations 
will yield, by virtue of the laws of large nunbers, a statistically 
meaningful description of the physical problem. 

Tte form of Monte Carlo used in this study is as follows: 


Define the system model by assuming 


a. 

b. 

c. 


d. 


model regression function 

method of error incorporation . , 

probability distributions of all errors (for all independent 

variables) _ . 

any equations used to model the phenomena of interest 


2. Use the computer and random sampling techniques to select 
values of the independent variables. 

3. Calculate dependent (output) variables using the prescribed 
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equations. 

4 . Estinate regression parameters for the assumed mdel. 

5. Replicate the experiment, each time with a new set of input 

values. 

Use appropriate statistical methods to calculate properties of 
the distribution of parameter estimates. 


6 
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E. Brief Sumrrary of Results 

A ply node fro. the AS-O-aphite /IHB epoxy oonposit. system is 
studie d. The nonte carlo scheae is used to generate a n»*er of 
response results, -rich are analyzed in graphical and nodical for* to 
supply a random process description of ooaposite ply elastic constats, 
theranl expansion coefficients, and strengths. Histogram and 
distribution plots of results for assured narrow and wide oru-iations in 
input properties are cospared with a deteradnistio base case for an 
aligied ply. The figures denonstrate the range of oal«s that re po 
variables assure for the exasple data under consideration. 

Confidence internals are calculated for response variables in 
subsequent sanples, which are nornelized with respect to an appropriate 
independent variable, to yield plots of normalized response as a 
function of fiber voluee ratio, for various val«s of distribution 
paraeeters for the related independent variable. These plots 
denonstrate the sensitivity of ply properties to random selected 
uncertainties in constituent and fabrication variables. 

several -iltipl. linear regression nod. Is were calculated for 
response variables. The relate correlation of predictor (independent) 
variables with response is studied for all output properties consid^. 
Varying levels of significance **re achieved in the regression 
equations, d* to tie differed in conplexity of response -tables. 
Elastic constants can be described adequately with staple regressor 

functions, and genially explain tat— ~ — •» «— * **“ 

variations about a — . The regression ecdels 


observed response 
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studied far strength, although achieving better reliability with higher 
order regressor functions, demonstrate such low signif ioanoe as to be 
practically useless for predict ive purposes. This is not an unexpected 
result, because of the coiqp lex nature of strength behavior in composite 


materials. 



CHAPTER II 


PETHODG OF CALCULATION 


A. Overall plan 

1. Input structure for ICAN 

The input data for a typical execution of the available ICAH 

program consists of (see Fig 5) 

1. leader card 

2. control cards 

3. ply data cards 

4. iioterial system cards 

5. load cards 

For repeated use of the ICAN program, input data files mist be 
created and used one at a tine. Each successive run of the nester 
program (of viiich ICAN is nade a subroutine) writes the input file from 
user-supplied paraneters aid calls ICAN. Hie ply data cards contain 
randomly generated fiber orientation angle vali*s. The naterial system 
cards contain randomly generated valtes for fiber and void volune 

rat ios . 


17 
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2. Constituent Property Variations 

Each w*»essive execution of ICflN uses a distinct set of rater ial 
properties for fiber and ratriK. The random nonber generation is 
perfomed with user-supplied parameters nfcich are stored in a separate 
file. The options of using either generated properties or using the 
values contained in the resident data bank are available. Any subset of 
the parameters described ray be generated or held constant with proper 
specification of the Boo leans vfiich control the input to the ICAN 
program, (see Figs. 6,7) 4 
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FIBER STRENGTH VARIES 
STDATA 15 

T 50 T 

F 000.0 

F 
F 
T 

PLY 

MATCRDAS-IUIHS 
PLOAD 10. 

PLOAD 0.0 . 

PLOAD 0 . 0 

OPTION 0 


; CONSTANT FIBER VOLUME RATIO OF 0.30;TAPE 003131 


1 15 T 

F T T 

3.00 



10.0 0.300 0.200 

5 


70.00 70.00 

.0 

.57 

.000 

AS-1IMHS 

0.0 0.0 0.0 

0.0 o.o 

0.0 

.03 


0.0 


Fig. 6- Command Input 


EFP1 

EFP2 

GFP12 

GFP23 

SFPT 

SFPC 

EMP 

SKPT 

snpc 

snps 


0.3100E 08 
0.2000E 07 
0.2000E 07 
0.1000E 07 
0.6C00E 06 
0.4000E 06 
0.5000E 06 
0.1500E 05 
0.3500E 05 
0.1300E 05 


0.3000E 07 
0.2000E 06 
0.2000E 06 
0.1000E 06 
0.1000E 02 
0.1000E 02 
0.5000E 05 
0.1000E 02 
0.1000E 02 
0.1000E 02 


Fig. 7- Constituent Variation Input. Example for AS;} Graphite 
9 fiber and IMHS Epoxy matrix, with wide variations of 
stiffnesses and strengths. 
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3. Repeated runs 

The user must specify the number of I CAN runs desired in a given 
sample. In this study, fifty (50) runs were used throughout, to take 
advantage of the simplification in statistics by using suitably large 
samples. From elementary statistics, it is known that any process that 
is the result of the combined interaction of several probabilities can 
be assuned to approximate a normal distribution. For phenomena that are 
assumed to approximate a normal distribution, the simplest forms for 
calculating statistics apply to suitably large samples (usually greater 
than thirty). The sample size of fifty was chosen to supply a 
practicably large amount of data, within the restrictions imposed on 
computation time. 

The data generated by repeated execution of the ICAN routines is 
stored in a sequential access dataset, where the 5® output files are 
separated by end of file markers. This arrangement allows a single 
Fortran unit to be used for output throughout. A simple fiovchart of 
the data generation routines is shown in Fig. 3(a). 

4. Data collection 

The ICAN output files are searched to locate the specific material 
properties and strengths of interest in this study. The flovchart of 
data collection routines is shown in Fig. 8(b). After obtaining the 
sample of ICAN output, the investigator may choose to scrutinize 
parameters or calculate statistics aside from those chosen in this 
study. This is likely, in light of the large quantity of data available 
and the t«ed for limiting the scope of this particular study to 
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representative properties. The user sould have to s^ply additional 
cede or adapt eaisting node to suit his purposes in this oase. The 
coded codifications to ICftK used in this sttrfy are inclaied in dppendis 


A. 
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B. Generation of Pseudo Random Nunbers 

An integral part of any rente carlo simlation la tie ose of rardo. 
carters having a spmified distribution which is assured to characterize 
the process under study. Indeed, reny statistics terthooks carry tables 
Of random nuahers as appendices. Sirulations using large sarples 
require many repeated calculations, each with different "random" 
nnrters. Since filling of a ccputer remry with a large table of 
random nunbers is wasteful, algorithm have been developed («ef. 9) to 
generate streams of random nunbers whenever needed in the process of 
calculations. The nurhers used are usually obtained using sore form of 
a recursion relation, hence the sequence is terred pse«Io-random. 


1. Uniform Distribution 

The starting point for rany random nunfcer scheres is the uniform 
random neater generator, which sirulate. a sarple from the uniform 
distribution. A continw«»us random wmriatle has a miferm distribution 
«er an interval . to b ( h > a ) if it i. *sually likely to t^we on any 
valwrn in this interval. The probability density fimetien is thus 
constant over ( a,b ) and has the form 


f(-) * -b^T * 5 * 2 “ 

= ® elseniiere 


The probability distribution function is, on integrating 

F(x) = 0 x { a 


< x < b 


b - a 
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= 1 x > b 

The uniform distribution is shown in density and distribution form in 
Figs. 9a and 9b. 

Lefaner (Ref. 1®) proposed the congruent ial aethod of generating 
p seudo random numbers conforming to the uniform distribution. The 
recurrence relation takes the form: 

x. = (ax. . + b ) modulo m 
x v l-l ' 

vtere the notation signifies that x A is the reminder viien (ax^ * b ) 
is divided by m. The multiplier a, increment b, and modulus m are 
integers. The starting value x 0 nust be assuned, and is knovn as the 
"seed" of the generator. Generators for vhich b = 0 are knovn as 
multiplicative. Hey are called mixed vhen b is nonzero, fecause 
selection of tte multiplier a and modulus m strongly influence the 
generator, most generators in use are of the multiplicative form. A 
discussion of the choice of parameters, maximum period, and degree of 
correlation of this generator is available (Ref. 11). 

For a given uniform random number u on the interval (0,1) a random 
number x having a desired distribution F(x) is often obtained by solving 
the eqmtion u = F(x) for x (Ref. 12). Since the process requires the 
determination of the inverse distribution function F *(x), its use 
depends on the ease of deriving the expression or some approximation. 

The following sections describe the distributions used, and methods for 
gemrating random numbers on those distributions. 
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2. lfariiHl (Gaussian) Distribution 

The Host oonrnm distribution is the familiar nornwl distribution, 
with the "bell shaped” density function, given by 

1 f (x-M) 2 1 

-<» < X < co, fi < «, and a > 0 

with nean fi and standard deviation a. The distribution function is 
written 


1 

X 

(U-/0 2 

F(x) * JSS* 

exp 

-CO 

■ ~1P~ 


which cannot be expressed in closed form analytically but can be 

numerically evaluated at any value of x. 

■n* Box-Miller or "Polar" netted (Bef . 13) is test coteonly used 
for generating ration, deviates from a tean to approxiwte the nortel 
distribution. If Xj and x., are independent uniform random variables, 

then 

y t = a [-2 In x 1 )®' 5 cos 2nx 2 + fi 
y 2 = o(-2 In x/- 5 sin 2nx 2 ♦ fi 

are independent random variables with the standard nornel distribution 
aixi staixiard deviation a. 


having nean fi 
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3. Ganna Distribution 

U* gansia distribution is a two-parameter distribution is 

flexible in fitting a variety of random processes. It is a one sided 
distribution in that physical quantities that are limited to valtes in 
the positive range are freqiently nodeled by it. Its density function 

is giwn by 


f(x) = 


-Xx k-1 
? x 


r(k) 


Were x, X, k > 0, and k is an integer. 

The paraneters X and k nay be interpreted as scale and shape paraneters, 

respectively. r( k) is the veil known gamra function, 

, Too k-1 -u . 

f( k ) - Jo u e du » 

vfcich is widely tabulated. The gatma distribution function is given by 


F(x) = 


T(k) 


% 

X 

k-1 -Xu 
_ u e 
0 


du 


r(k,xx) 


x > 0 


= r(k) 

s 0 elsexdiere 

vtere r(k,u) is the inoonplete gamna function 

T(k 


* 

x U k_1 
t,u) = x 

0 

* 


-K. 

e dx 


vdiich is also widely tabulated. For integer values of k, 

T(k) = (k-1)! 

and the ganaa distribution is knowi as the EWangian distribution after 


A. K. Erlang, %*o introdt«ed it in the theory of que«s and Ifarkov 
processes. 
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Ganna variates are generated using the seqiKice 

U 1 ,U 2 ,U 3* °k 

satisfying the uniform distribution on the interval (0,1). 


The recursion relation is 


- - 


— ln V 




vfiere x is a gamm variate having parameters X and k (Ref. 14) 
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4. Weibull Distribut ion 

The Vfeibull distribution (Ref. 15) is Host popular **e n node ling 
problem of reliability, imterial strength, and fatigi*. The Ifeibull 

density function is qiwn by 

f (x;<x,0) = a^ _1 exp(-uit^) 

® < x < «, a>®, P>1 

where a and fi are the shape and scale parameters, respectively. The 

cumulative distribution function 

y = F(x) = 1 - exp[-(K/£) a ] 

leads iimediately to the inverse relationship 

F -1 (y) = x = - £[ ln(l-y) ] 1/a 

as tto desired Weibull random generator vhen y is a uniform random 
variable. 

Figures 9-12 stow tte above distributions in analytical form. 
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q Distribution Assunptions 

^ variables chosen tor variation are those tor -.ich reason*.!, 
assertions can he »de to describe their distribetion. The tiber 
geonetric configuration with respect to ply aaes is assmd to folic- a 
^1 distribution snth sean of zero (degrees) and sor saell standard 
deviation, to be specified. The fiber volune ratio is assesed to be 
norselly distributed about sone seen between 8.3 and 8.7. The void 
eelur ratio, w*i=h is ideally ssell, is assurd to follow a gam 
distribution skewed toward zero. (Note that in the P- distribution 
used, a unite of zero has a probability of zero. This sodel is chosen 
begause tie state of rs, present tenufactieinp technolo*. precluies the 
fabrication of a fiber coeposite cospletely free of void.) 

The properties of irdividual fibers and setri. are varied. The 
uueaal and shear aoduli are assueed to fcilow the norsal distribution, 

and the strengths are assused to be Ueihull distributed. 

Figs. 13-27 show the results of randon nesber generation in each 

distribution studied. The density (or histo^as.) and curative 
distribution plots are show*. Several ^eibull and ganro distri 
sintilat ions are sKm. to desonstrate the effects of assmd paraseter 
the distribution sanpling. 


variations on 
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HISTOGRAM FOR 

normal generator 



distribution OF 
NORMAL GENERATOR 



Fig. 


13- Normal Distribution Simulation with mean of 
0.0 and standard deviation of l.u* 
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histogram FOR 
GAMMA GENERATOR 



distribution OF 
6AMHA GENERATOR 
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(a) histogram 


distribution OF 
GAMMA GENERATOR 
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DISTRIBUTION 

HEIBULL GENE! 


is to gran 
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histogram FOR 
WE1BULL GENERATOR 



distribution OF 
WEIBULL GENERATOR 
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histogram for 
heibull generator 



distribution of 
heibull generator 
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DISTRIBUTION OF 
HE I BULL GENERATOR 




histogram for 
he I bull generator 
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distribution of 
heibull GENERATOR 
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histogram for 
heibull GENERATOR 



distribution of 
heibull generator 




47 


D. Use of ICflH 

This section describe, the essential theories — .scions 
imrporated in «. ** P^- »• 

(emulations, a* definitions are inclined in «pp«s.ix B. 


1. Conposite Hicronechanics 

the branch of opposite seohanios -uoh relates ply properties to 
constituent properties is as -resit. ^ioroseohanics. H- inputs 

consist rot only of constant -teri.l Properties (fiber and 
but geonetric configuration - — tion P™ 

ply hygral , therm 1, ard eechanical pmperties. Tbe assertions for 

equation development are: (Ref. 16) 

The Mechanics of feter identlf 

^Tl^mtHKiLe loads according to the sohemtio 

^ifaS’m^stituents behaue in a U— 

^e f p^rt^s^-lf isotropic in the 2-3 plane. 

Jc£ eaists^at'the fibermetri. interface. 


1 . 

2 . 

3. 

4 . 

5. 

6 . 


. . ^ terminology used in the equations 

The direction conventions and term ogy 


ares 


1 . 

2 . 

3. 

4 . 


Properties seasoned along fiber direction are callad 
SS^erieasured transverse to fiber direction are called 

S^siear is a lsc tf^y .trial 
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2. Laminate Theory 

Classical laminate theory supplies a convenient prooedtre to 
predict the response of a laminate to external load. The theory uses 
anisotropic elasticity to obtain the stress-strain relationship for the 
basic lamina. The stress-strain relations of individual laminae are 
transformed to coincide with a global set of reference axes. The 
stress-stain law of the laminate in terms of the properties and 
distribution of individual laminae are calculated using a sunnation. 
Resultant forces and xonents are defined by integrating the stresses 
through the thickness of the laminate. The plate constitutive equation 
is inverted, giving midplane strains and plate curvatures in terms of 
applied forces and nonents. These strains and curvatures are 
substituted into the lamina stress-strain equation to obtain lamina 
stresses in the global system. The stresses obtained are then 
transformed into the principal material system of the lamina in question 
and conpared with ultimate stresses obtained using failure criteria. 
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3. Strength Theories 

• Tran sake use of several assunptions. 

•He strength theories in ICflH *eke use o 

, that there are five characteristic val«s of 

First, it is assuned that there 

strength of a unidirectional conposite: 

1 . longitudinal tensile strength 

2 . longitaiinal co.p.essioe strath (3 separate criteria) 

a. rule of mixtures 

b. fiber microbuck ling 

c. de lamination 

3. transverse tensile strength 

4. transverse conpressive strength 

5 in-plane or intralaminar shear strength 
^ fracture nodes usually associated with these strati* are sho~ 

schenat ical ly in Fig. 29. 

ply strengths are caicuiated (in the pi, coordinate syste-). 

geonetric transitions are used to oaiouiate conposite — ioads. 

Tte process us«I is briefly described below. 

, calculate lc-d. (in oolite system) reguir- to ind« load egua 

to ply strengths (in ply systems) for each .ode. 

2 . calculate -ini™, of failure to*. 

3. calculate minium of failure loads of all plies, *™1 «« th,S 1 
„ the failure strength of the ctmposite for a p^icular failure 


node. 



STRAW 


Fig. 28- 


Typical Stress-Strain behavior of 
unidirectional fiber composites. 



•1 unfRutftMl 
tonUon. 



mcrmucuinc 

ft) lanf RuftatJ camKHtlon. 



Flo 29- In-plane fracture modes of unidirectional 
S (ply) f1ber composites. 
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E. Heuiew of Applicable Statistical Crcepts 

Coapositc p~pect.es ace calculate* for — samples usrn, a 
specific set of distributions of iupot Pities. In this context, 

sx.il sanpling ««ry dees not apply, 
sufficiently large. 


1. Sanple ffeans 

Calculation of the *ean saxple ual*s proceeds by defining 


n 

L Ki 

Ttean = x — ji 

vfcere n = sample size 

x .= sanple data values 

* * «, «o tte sanple nean is assumed to be the 

The population nean is unknown, 

best estimator of the population nean. 

2. Sanple Standard Deviation 

estinste nf the population sta^rd delation is calculat 
using the statistically efficient estinntor 




-. 1/2 

, - *)* I , n > 3® 


3 Confidence Interval Estimates 

_ , in ^ area of statistical inference is the 
An important problem m the 

ters (st*h as nean, variance, etc.) 
estimation of population parameters 1* 

_ 3 and o are the nean and standard 

sample statistics. Parameters x and 
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deviation of the sanpling distribution of a statistic S. The saapling 

distribution of S is a* sound as approxixetely nonial (sfiioh is trs* for 

statistical distribution, if n 1 38). Confiifanoe interval 

estimtes are constructed for the statistic S. Thus, intervals as 

identified for shich it can be asserts with a reasonable de^ee of 

certainty that they contain the paraseter considered. Obviously, the 

degree of certainty (or confidence level) will vary with the size of the 

interval chosen. Values of confidence coefficients, * c , are associated 

with confidence levels. For exaeple, •» act®' saeple. statistic S is 

expected to be found lying in the interval (5 - to (x * z c o) (— " 

„ is the unknow. population standard deviation) we percent of the 

tixe. Let the z c val® in this exasple be 1. Bssunun, a "cruel 

saxpling distribution, (with z o = 1) the noreel distribution area 

fraction specifies that S falls betveen (5 - o) and (x * o) about 

£8.27/. of the tixe. Sinilarly, the confidence of x lying in the 

interval (S - o) to (S . o) is about 68.27/.. The endpoint, of the 

intervals are know. a. confidence lieits. Various confidence 

coefficients z . corre^onding to freqwmtly used oonfittenc. levels, 
c 

have b ee n tabulated. 

In this work, the confidence interval for neans is given in tern* 


of the sanple statistics by 


k + z 
- c 


raere z is the confidence coefficient, af.ich takes on values of 
1.645, 1^960, ed 2.588 for the 98, 95. and 997. confidww. levels. 
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respectively. 

4. Regression 

The term "recession" as used in the are. of statistics raters to 
the process ot for-.!. tic, a sethesatical sodel to eaplaio raodeeiy 
observed pta-sona. Sose factional tore tor the -y each oariahle 
enters the sodel sust be assumed. Comparison ot the degree ot tit ot 
ditterent assused sodels ideally leads to a better -del- *■ to5io 

regression strategy used here consists ot! 

1 . Assu* a eultiple li^ar regression sodel. Ihe nor*, actions 

for such a nodel are: 

m = rxio> + {*} 


where 

(Y) = vector of dependent variable vaix*s 

[X ] = matrix of functions of independent variable 

{£} = regression “trte” valtes 

{e} = errors 


The raariral equat ions can 


be solved as follows: 


vhere 


[x] t {y> - [x] T m<» * m t w 

<b) = [X T xf 1 [xT{V) 

^b} = parameter est ilia tes 


2. Use a standard statistical package 

paraneters • 


(Ref. 17) to estimte regression 
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t 

3. Calculate properties of regression parameter distributions to 

assess irodel precision. 

In the event that [X T X] is singular, implying that some of the 
nxmal equations are linearly dependent, [X T X] 1 does not exist. The 
nodel should be expressed in terms of fever parameters, or should 
irclude assured restrictions on the parameters. 

2 

Tte square of tie multiple correlation coefficient, R , is usually 
calculated for each regression model, and supplies a convenient measure 

VS 

of the degree of fit betveen data vali»s {Y} and valies {Y} 
predicted by the regression equation. It is defined by 

2 Sum of Squares die to regression model 
R = Total Sum of sqimres about mean Y 

T (Y. - Y) 2 
r (Y. - y)* 

Frequently, it is necessary to determine if inclusion of particular 
terms in a regression model is warthtfiile. To this end, the extra 
portion of tte regression sum of squares vhich arises dte to the terms 
under consideration is calculated. The mean square (defined as the sum 
of squares divided by the corresponding degrees of freedom) derived from 
this extra sum of squares can be compared with s s , the estimate of a*, 
to see if it appears significantly large. If it does, the terms under 
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consideration should be incited. The statistic is frequently compared 
to tie appropriate percentage point of the F- distribution, %faich is 

tabulated. 

Supopose the extra sum of squares die to a paraneter, given that a 
number of otter parameters are already in the model, is calculated. 

Syntxi 1 ica 1 ly , 

SS ( b i Jb 0 ,b l , “* ,b i-l ,b i+l’*' ,,b k ) * = 1,2 k 

represents a one degree of freedom ( 1 df ) sum of squares viuch 

treasures tie portion of the regression sum of sqieres due to the 

coefficient b^ This is a measure of the valie of adding a term to 

the nodel OUch previously did not include /T. The corresponding nean 

square, equal to the SS (siree it has one df) can be compared by an 

F- test to s 2 . This is known as a partial F- test for the single 

paraneter 0., which is a special case of the F- test described earlier. 

The stepwise regression procedure (Ref . 18) is a structured vay to 
insert variables in order of correlation until the regression equation 
is satisfactory. The partial correlation coefficient measures the 
relative importance of terms r»t yet in the nodel, to choose the next 
candidate for entry. The analagous statistic, F- to enter (or F- to 
remove) is usually evaluated for each predictor at every stage as though 
it were the last term to enter the model, to determine if terms retained 
at a previous step have become superfluous, because of some linear 
dependence with terms now in the nodel. The largest F- statistic 
calculated at each step is compared with the appropriate percentage 
point of the F- distribution, and the predictor variable is entered (or 
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renewed) based on the signif icance o£ this F- test. Testing of the 
least useful predictor is perfomed at every step. The R* statistic is 
calculated, to provide a measure of the val« of the regression at each 
' step. This stepwise linear regression scheiie is used in this «3rk 
because of its computational ecanony, and because it allova the analyst 
to assess the relative influence (or correlation) betveen individual 
predictor variables of a selected nodel and response for a particular 
data sanple. Otter schemes are available (Ref. 18 ) » such as backward 
elimination. The stepwise procedure is recoimended for its direct 
nature in testing the nodel with only significant predictor terms. 



CHAPTER III 


RESULTS 


A. Property Histograms and Distributions 

In this snr it, fiber and untris properties are aliased to assone a 
range of valimrs to assess the sensitivity of the opposite ply 
properties to constant perturbations, a-aphite fiber and epoxy 
aatrix are used as tte constituents. Initially, t«. separate samples of 
output data are generated and studied to demonstrate the effects of 
input paraneter charges on composite material properties. These tvo 
oases are compared with a deterministic base case with no random input 
property generation. The data for all three oases is given in Table I. 

The results of cases 2 and 3 are sbovar in histogram and cumulative 
distribution form in Figs. 30 - 42. The resuit, of the deterministic 
oase 1 are summu-ioed in Table II, and can be easily compared with the 
histograms and distributions. 
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tqctf T- IMPUT DATA FOR SfifTLIHG 


INPUT 

CASE 1 

CASE 2 

CASE3 

THETA (degrees) 

©.0 

- 

- 

fi 

— 

©.© 

0.0 

a 

- 

5.0 

10.0 

FVR 

©.5© 

- 

- 

fi 


0.5 

0.5 

a 

- 

0.1 

0.2 

WR 

©.©1 

- 

— 

X 

- 

3.0 

3.0 

k 

- 

3 

5 

EFPl(ksi) 

3100© 

— 

— 

fi 

- 

31000 

31000 

a 

- 

1500 

3000 

EFP2(ksi) 

200© 

— 

— 

fi 

- 

2000 

• 2000 

a 

“ 

100 

200 

cari2(ksi) 

200© 

— 

— 

fi 

- 

2000 

2000 

a 

- 

100 

200 

GFP23(ksi) 

10M 

— 

— 

fi 

- 

1000 

1000 

a 

- 

50 

100 

SFPT(ksi) 

4G0 

- 

— 

fi 

— 

400 

400 

a 

- 

20 

10 

SFFC(ksi) 

400 

- 

— 

fi 

- 

400 

400 

a 

- 

20 

10 

EJF(ksi) 

500 

- 

— 

fi 

• 

500 

500 

a 

- 

25 

50 

STPT(ksi) 

15 

— 

— 

fi 

— 

15 

IS 

a 

- 

20 

10 

SfPC(ksi) 

35 

— 

— 

fi 

— 

35 

35 

a 

- 

20 

10 

SfFS(ksi) 

13 

- 

— 

fi 

- 

13 

13 

a 

- 

20 

10 


toct F ii- CASE 1 RESULTS 


PROPEHIY 


VALUE 


ECU 

EC22 

EC 12 

KUC12 

KUC21 

CTE11 

CTE22 

CTE12 

SCXXT 

SCXXC 

SCYVT 

SCYYC 

SCXYS 


15750 ksi 
1065 ksi 
516 ksi 
0.275 
0.018 

0.775 x 10~ J 

0.181 x lO 
0.000 
203 ksi 
165 ksi 
11.74 ksi 
27.41 ksi 
10.01 ksi 


CUMULATIVE FftEOUCMCV 







CUMULATIVE 


HISTOGRAM FOR EC22 

transverse modulus 


I DU • 1 

RANGE t c 07> 


(a) case 2 histogram 


histogram for 
transverse modulus 


(b) case 3 histogram 


DISTRIBUTION OF EC22 
T8ANSVEBSE KOOULUS 


DISTRIBUTION OF EC22 
transverse modulus 


IHJ 


RANGE < F 07> 

(c) case 2 distribution 


ftAM&E C E 07> 

(d) case 3 distribution 


^C; case l — 

Sampling results for Transverse Elastic Hodulus 
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HISTOGRAM FOR EC12 HISTOGRAM FOR EC12 


SHEAR HOOUIUS 


SHEAR HOOUIUS 



DISTRIBUTION OF EC12 DISTRIBUTION OF EC12 

SHEAR HOOUIUS SHEAR HOOUIUS 



Fig. 32- Sampling results for In-plane Shear Modulus 


CUMULATIVE 



RANGE 


< E 00> 


RANGE 


(a) case 2 histogram 


(b) case 3 histogram 


DISTRIBUTION OF NUC12 
POISSON RATIO- MAJOR 


DISTRIBUTION of NUC12 
POISSON RATIO- MAJOR 




RANGE < E 00> 

(c) case 2 distribution 


RANGE < E 005 

(d) case 3 distribution 


Flo. 33- Sampling results for Poisson Ratio (major) 




CUMULATIVE FREQUENCY 





CUMULATIVE FBEOUCMCT 






CUMULATIVE 
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HISTOGRAM FOR CTE22 
TRAN. THERMAL EXPANSION 
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HISTOGRAM FOR CTE22 

TRAN. THERMAL EXPANSION 



(a) case 2 histogram 


(b) case 3 histogram 


DISTRIBUTION OF CTE22 DISTRIBUTION OF CTE22 

TRAN. THERMAL EXPANSION TRAN. THERMAL EXPANSION 



Fig. 36- Sampling results for Transverse Thermal Expansion 
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HISTOGRAM FOR CTE12 HISTOGRAM FOR CTE12 



DISTRIBUTION OF CTE12 DISTRIBUTION OF CTE12 



Fig. 37- Sampling results for Thermal Expansion Coupling 
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histogram for scxxt 

LONG. TENSILE STREN6TN 


range 


( E 03> 


(a) case 2 histogram 


histogram FOR SC) 

LONG. TENSILE STRENGTH 



5 ?. 5 


RANGE 


< E 03) 


(b) case 3 histogram 


distribution of scxxt 

LONG. TENSILE STRENGTH 


DISTRIBUTION OF SCXXT 
LONG. TENSILE STRENGTH 




RANGE < t 03) **»■« * 1 03 > 

(c) case 2 distribution (<J) case 3 distribution 

a 38- Sampling results for Longitudinal Tensile Strength 
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histogram for scxxc histogram for scxxc 



distribution of scxxc 


DISTRIBUTION OF SCXXC 



Fig. 39- Sampling results for Longitudinal Compressive Strength 


CUMULATIVE FREQUENCY 


70 


HISTOGRAM FOR SCYYT HISTOGRAM FOR SCYYT 

TRAN* TENSILE STRENGTH TRAN. TENSILE STRENGTH 



RANGE < E 02) MAM6E t E 02> 

(a) case 2 histogram (b) case 3 histogram 


DISTRIBUTION OF SCYYT DISTRIBUTION OF SCYYT 

TRAN. TENSILE STRENGTH TRAN. TENSILE STRENGTH 



RANGE <.E 02) RANGE < £ 02) 


(c) case 2 distribution (d) case 3 distribution 

Fig. 40- Sampling results for Transverse Tensile Strength 
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histogram for scyyc 
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histogram for scyyc 
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(a) case 2 histogram 


RANGE 
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(b) case 3 histogram 


1STR1BUT10N OF SCYYC 
TRAN. COMPRESS STRENGTH 


distribution of scyyc 

TRAN. COMPRESS STRENGTH 


RANGE < t 02> 

(c) case 2 distribution 


RANGE < E 02) 

(d) case 3 distribution 


41- Sampling results for Transverse Compressive Strength 
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HISTOGRAM FOR SCXYS HISTOGRAM FOR SCXYS 

IN-PLANE SHEAR STRENGTH IN-PLANE SHEAR STRENGTH 



(a) case 2 histogram (b) case 3 histogram 


DISTRIBUTION OF SCXYS DISTRIBUTION OF SCXYS 



Fig. 42- Sampling results for In-plane Shear Strength 
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B. Fiber Strength Effect 

To show tte effect of fiber strength changes on the langitxxiinal 
strengths of the conposite, several shape paraneters of the weibull 
iistribution for fiber strength are assumed. The iionte carlo procedure 
is then conducted at several fiber volune ratio values. All properties 
are varied, except fiber volune ratio. The distribution paraneters of 
all properties except fiber strengths are held constant. The curves 
generated are shove in Figs. 43 and 44. In the figures the solid lines 
and syntools show the neans of the 957. confidence interval estinetes for 
the sanple size of 5© chosen at each point. The points on both sides of 
each curve locate the upper and lover bounds of the confidence 
intervals. The convention described is intended to provide a convenient 
indication of the dispersion of the sanple valies at each point. 



NORMALIZE!) (SCXXT/SFPT) < E -02) 


74 


LONG. TENSILE STRENGTH 



Fig. 43- Longitudinal Tensile Strength; for various 
shape parameters of fiber strength. 


NORMALIZED (SCXXC/SFPC) < E -OE) 


LONG. COMPRESS. STRENGTH 



Fig 44- Longitudinal Compressive Strength; for various 
S ’ shape parameters of fiber strength. 
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C. ttrtrix Strength Effect 

Tl* effects of changes in netrix strength on oonposite strengths 
are studied by suitable variation of the shape parueters gowrning the 
■etrix strength distributions. Analagous to the plots gi\«n for fiber 
strength effects, the matrix effects are showi in Figs. 45 - 47. 



NORMALIZED <SCYYT/SMPT> < E -02) 
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TRANS* TENSILE STRENGTH 



FIBER VOLUME RATIO 




NORMALIZE!) (SCYYC/SMPC) ( E -02) 
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TRANS. COMPR. STRENGTH 



FIBER VOLUME RATIO 

Fig. 46- Transverse Compressive Strength; for various 
shape parameters of matrix strengths. 



(SdUS/SAXOS) QaznVHttQN 
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IN-PLANE SHEAR STRENGTH 



Fi«. 47- In-plane Shear Strength; for various 
3 shape parameters of matrix strengths. 
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D. Fiber Orientation Effect 

Assum'd val« of tie fiber orientation angle distribution 
paraneter are consecutively used in the sente carlo procedure to assess 
the effects on several oonposite properties. These plots are shovn in 

Figs. 48 - 57. 

E. Fiber Stiffness Effect 

Assuned values of the fiber modulus distribution paraneter are used 
in tl* simulation to similarly assess the effects on the related 
conposite properties. The plots thus generated are shovn in Figs. 


58-67. 



NORMAL 1 ZED <ECII/EFP1> 
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10.0 
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Fig. 48- 


Longitudinal Elastic Modulus; for various 
shape parameters of fiber orientation. 








NORMAL I ZEO (SCXXT/SFPT) < E -02) 
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LONG. TENSILE STRENGTH 



FIBER VOLUME RATIO 

Fiq. 51- Longitudinal Tensile Strength; for various 
shape parameters of fiber orientation. 



NORMAL I ZED <SCXXC/SfrPC> < E -02) 




(20' 1 > UdMS/UAOS) 03ZnVHtt0N 
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TRANS. TENSILE STRENGTH 



FIBER VOLUME RATIO 

Fig. 53- Transverse Tensile Strength; for various 
shape parameters of fiber orientation. 



NORMAL 1 ZED <SCT YC/SHPC) < E -02) 
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TRANS. COMP. STRENGTH 



NORMAL I ZEO (SCXYS/SMPS) 
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POISSON’S RATIO (MAJOR) 



Fla 56- Poisson’s Ratio (major); *>r various 
9 * shape parameters of fiber orientation. 





LONG. ELASTIC MODULUS 



Fig 58- Longitudinal Elastic Modulus; for various 
9 shape parameters of fiber modulus. 



TRANS. ELASTIC MODULUS 



fiber volume RATIO 

iq. 59. Transverse Elastic Modulus; for various 
shape parameters of fiber modulus. 


IN PLANE SHEAR MODULUS 



fiber volume ratio 

. 60- In Plane Shear Modulus; for various 

shape parameters of fiber modulus. 
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POISSON’S RATIO (MAJOR) 



FIBER VOLUME RATIO 

Fig. 61- Poisson's Ratio (major); for various 
shape parameters of fiber modulus. 


NUC2 1 < E -03) 
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POISSON’S RATIO (MINOR) 



Fig. 62- Poisson's Ratio (minor) for various 
shape parameters of fiber modulus. 


uo- 3 > (idjs/ixxos) miiummom 


LONG 


TENSILE STRENGTH 



F< _ 63- Longitudinal Tensile Strength; for various 

Fl9 ‘ shSpe parameters of fiber modulus. 
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TRANS. TENSILE STRENGTH 
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FIBER VOLUME RATIO 

Fin. 65- Transverse Tensile Strength; for various 
shape parameters of fiber modulus. 



NORMALIZED (SCXYS/SMPS) < E -on 
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IN PLANE SHEAR STRENGTH 



Fin. 67- In Plane Shear Strength; for various 
shape parameters of fiber modulus. 
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G. Regression ltodels 

The CHrtpat data =t cases 2 thresh II are T«d a. 5*oes.iw 
to the regression schere. The goal of stepwise regression, as osed 
tow, is to Treasure the degree of correlation betreen a dependent and a 
set of independent oariables for a gioen set of data. Hie outputs of 
the regressions eondreted show the independent variables accepted into 
the redel (based on F-test criteria) in order of degree of correlation 
with the dependent variable of interest, along with the final R 
statistic, (the R 2 valres represent the square of thereiltipl. 
correlation coefficient, a convenient Treasure of the fit between data 
vaiues and valres predicted by the regression equation.) 

The ordering of predictor variables by stepwise regression has 
several inportant uses. In this sttrfy, the schese facilitates easy 
investigation of the effects of eeterial changes on conposite 
properties. Since the rente carlo schese permits generation of large 
areunts of data, the regression is easy, inespensive, and ere provide 
insight concerning the sensitivity of dependent vreiable. for as.»ed 
distribution, of predictor variables. * variety of aeterial 
^figurations and constituent distributions are esasaned, and a redel 
constructed for each depe«len« (or response) vreiable. It rest be noted 
that the relative correlations of predictor variables with response 
variables will be functions of tie assured distributions, the preticular 

v s r functional ranner in tfnch the predictor 

data sanple considered, and the tuncx ion- 

variables are incorporated into the wodel. 

» siaple regression redel was assured for each response variable. 
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The first set of "simple" regression node Is uses as predictor functions 
only the iwlependent variables as individual terms. To be sore precise, 
the predictor variables used are not simply the independent variable 
values, for there are 15 of these for each layup. Hie arithmetic scan 
of independent variable values is thus used as the predictor variable in 
the first set of regression models. Hie only exception to this is the 
use of the sin* of the average of the fiber orientation angles as the 
angular dependence predictor, denoted by THETft in the tables to follow. 
The simpler response variables can be adequately described using the 
linear function forms in the regression nodels. Hie simple variables 
include the elastic constants, (BC11, EC22, EC12, NUC12, NUC21) and 
coefficients of thernel expansion (CTE11, CTC22). Hie results of the 
regressions performed in the “simple" manner are given in Tables III - 
XIV. In the tables the input labeled with HI throtxjh H5 and W1 through 
U5 represent narrow and wide distributions of all properties. Input 
labeled N6 through H10 and U6 through U10 describe the same 
distributions, except that the composite is assumed unidirectional, x.e. 
no angular variation. Hie distinction shows the rediction in predictive 
capability induced by deviations of the fibers from aligned orientation. 

The models assumed for the response (output) variables are of the 

form 

Y = B 0 + BiXj ♦ BjX 2 + B 3 X 3 + ... + B^ 

Y s response variable (ECU, EC22, EC12, etc.) 

B = regression parameters to be obtained 


viiere 



103 


X q = average of independent variable values through the 
thickness of the ply (THETA f FVR, WR, etc. ) 

Each nodel postulated contains all independent variables that 
appear in the equations for the related ply property (see Appendix B). 


v 



TQPf b* ttt— LONGITUDINAL M3C0US (ECU) 
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SITPLE MODEL 


INPUT 

FVR 

TESTE ACCEPTED 

HI 

0.3 

FVR.EFPl 

N2 

0.4 

FVR.EFPl, THETA 

N3 

0.5 

FVR.EFPl, THETA 

N4 

0.6 

FVR.EFPl .THETA 

H5 

0.7 

FVR.EFPl 

VI 

0.3 

FVR, THETA, EFP1 

U2 

0.4 

FVR.EFPl .THETA 

V3 

0.5 

FVR, THETA, EFP1 

W4 

0.6 

FVR.THETA.EFPl 

U5 

0.7 

FVR.EFPl, THETA 


R* 


83.17 

92.63 

94.02 

94.59 

84.00 

64.49 

89.88 

72.85 

65.37 

57.83 


N6 

0.3 

N7 

0.4 

N8 

0.5 

N9 

0.6 

H10 

0.7 

U6 

0.3 

U7 

0.4 

VB 

0.5 

V9 

0.6 

V10 

0.7 


FVR,EFP1,EIF 

FVR.EFPl 

FVR.EFPl 

FVR.EFPi.ETP 

FVR.EFPl 

FVR.EFPl.WR 

FVR.EFPl 

fvr.efpi 

FVR.EFPl 

FVR.EFPl 


99.83 

99.81 

99.69 

99.74 

99.77 

99.13 

98.40 

98.90 

99.59 

99.34 
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INPUT 

N1 

M2 

K3 

M4 

N5 

U1 

U2 

U3 

W4 

U5 

N6 

N7 

N8 

N9 

N10 

U6 

W7 

UB 

U9 

U10 


tqct g TV- TRANSVERSE FCOULUS (EC22) 


FVR 

SIMPLE IEDEL 


R 2 

TE3?TS flCUfcKlt-U 



0.3 

FVR, EFP 2 


83.50 

0.4 

FVR 


85.23 

0.5 

FVR, EFP 2 


91.83 

0.6 

FVR.EFP2 


93.26 

0.7 

FVR, EFP2, THETA 


93.06 

0.3 

FVR, THETA ,EFP2 


78.36 

0.4 

FVR, THETA, EFP2 

• 

90.73 

0.5 

FVR, THETA, EFP2 


80.15 

0.6 

FVR, THETA, EFP 2 


86.05 

0.7 

FVR, THETA, EFP 2 


87.14 

0.3 

FVR,EFP2 


87.13 

0.4 

FVR,EFP2 


86.15 

0.5 

FVR, EFP 2 


90.97 

0.6 

FVR, EFP 2 


93.47 

0.7 

FVR, EFP 2 


92.05 

0.3 

FVR, EFP 2 


79.72 

0.4 

FVR, EFP 2 


70.71 

0.5 

FVR, EFP 2 


81.92 

0.6 

FVR.EFP2 


88.62 

0.7 

FVR,EFP2 


84.05 



TftgLE U- SHEAR MODULUS (EC12) 


SIWLE MODEL 


INPUT 

PUR 

TERTE ACCEPTED 

R a 

N1 

0.3 

THETA, FUR, GW 

97.01 

N2 

0.4 

THETA, FUR, GW, GFF23 

98.85 

M3 

0.S 

THETA, FUR, GPP , GFP12 

97.50 

N4 

0.6 

THETA, FUR, GW 

98.01 

N5 

0.7 

THETA, FUR, GW , GFP12 

98.42 

U1 

0.3 

THETA, FUR, GW 

94.79 

U2 

0.4 

THETA, FUR 

94.27 

U3 

0.5 

THETA, FUR, GFP23 

93.71 

m 

0.6 

THETA, FUR 

95.62 

U5 

0.7 

THETA, FUR, GW, GFP23 

96.67 

N6 

0.3 

FUR, GW 

97.66 

N7 

0.4 

FUR,GW,GFP12 

98.02 

M8 

0.5 

FUR , GW , GFP23 

96.65 

H9 

0.6 

FUR,OP,GFP12 

97.11 

N10 

0.7 

FUR,GW,<FP12 

98.55 

U6 

0.3 

FUR,GW,GFP12 

96.93 

W7 

0.4 

FUR,GW,GFP12 

92.45 

UB 

0.5 

FUR,GW,GFP12 

95.16 

1)9 

0.6 

FUR, GW 

97.18 

W10 

0.7 

FUR,GW,GFP12 

96.90 



107 


T PHI- VI- POISSON’S RATIO, MAJOR (HUC12) 
SIJPLE MODEL 


INPUT 

FVR 


R 2 

N1 

0.3 

THETA, EFP1 

96.39 

N2 

0.4 

THETA, FVR 

97.88 

N3 

0.5 

THETA, FVR 

96.60 

N4 

0.6 

THETA, FVR, EFP1 

98.32 

N5 

0.7 

THETA, FVR, EFP2 

96.62 

U1 

0.3 

THETA, EFP1 

88.43 

U 2 

0.4 

THETA, FVR 

84.62 

W3 

0.5 

THETA 

89.48 

m 

0.6 

THETA, WR 

84.05 

U5 

0.7 

THETA, FVR 

92.05 

N6 

0.3 

FVR 

97.83 

N7 

0.4 

FVR 

98.48 

N8 

0.5 

FVR 

97.77 

N9 

0.6 

FVR 

98.40 

N10 

0.7 

FVR 

99.17 

U6 

0.3 

FVR 

97.32 

U7 

0.4 

FVR, WR 

96.45 

UB 

0.5 

FVR 

96.38 

U9 

0.6 

FVR,GFP12,EFP2 

98.34 

U10 

0.7 

FVR,EFP2 

96.96 
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TABLE VII- POISSON’S RATIO. ICHOR (MUC21) 


STIFLE IEDEL 


INPUT 

FUR 

TERIE ACCEPTED 

R J 

Ml 

0.3 

THETA, FVR 

91.15 

N2 

0.4 

THETA, FVR, EFP 1 

94.78 

M3 

0.5 

THETA, FUR 

94.31 

N4 

0.6 

THETA , FUR, EFP1 , EFF2 

97.18 

N5 

0.7 

THETA , FUR , EFP 1 

95.87 

U1 

0.3 

THETA, FUR 

90.87 

W2 

0.4 

THETA, FUR, EFP2 

89.86 

W3 

0.5 

THETA, FUR 

91.93 

W4 

0.6 

THETA 

92.57 

U5 

0.7 

THETA, FVR, EFP 1 

94.78 

N6 

0.3 

FUR,EFP1,EFP2 

95.64 

H7 

0.4 

FVR,EFP1,EFP2 

94.90 

N8 

0.5 

FVRjEFPl ,EFP2 

95.40 

N9 

0.6 

FVR,EFP1,EFP2 

93.12 

N10 

0.7 

FVRjEFPl ,EFP2 

91.83 

U6 

0.3 

FVRjEFPl , GFP12 

87.73 

W7 

0.4 

FVRjEFPl ,EFP2 

85.06 

U8 

0.5 

FVR,EFP1,EFP2 

84.29 

W9 

0.6 

FVR,EFP1,EFP2 

90.37 

W10 

0.7 

EFP1 ,FVR,EFP2 

91.42 
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INPUT 

N1 
N2 
H3 
N4 
N5 
W1 
U 2 
U3 
U4 
U5 

N6 

H7 

K8 

N9 

N10 

U6 

W7 

U8 

W9 

W10 


TQPfF UTII- LONG. THERM. EXPANSION (CTE11) 
SIFFLE rDDEL 


FVR 

TERTE ACCEPTED 

R 2 

0.3 

FVR, THETA, EFP1 

90.29 

0.4 

THETA, FVR.EFPl.WR 

94.46 

0.5 

FVR, THETA, EFP1 , WR 

95.72 

0.6 

FVR, THETA, EFPl.WR 

95.23 

0.7 

THETA, FVR 

87.63 

0.3 

THETA, FVR, EFP1 

80.53 

0.4 

THETA, FVR 

78.91 

0.5 

THETA, FVR 

84.77 

0.6 

THETA, FVR, WR 

74.37 

0.7 

THETA 

80.50 

0.3 

FVR,EFP1 ,WR 

97.21 

0.4 

FVR,EFP1 , WR 

96.96 

0.5 

FVR.EFPl.WR 

96.53 

0.6 

FVR.EFPl 

96.60 

0.7 

FVR.EFPl 

96.24 

0.3 

FVR.EFPl 

91.60 

0.4 

FVR.EFPl 

90.88 

0.5 

FVR.EFPl 

91.55 

0.6 

FVR.EFPl.WR 

96.03 

0.7 

FVR.EFPl 

94.13 



Tpmr IX- TRAMS. THERM. BjPAjSICM (CTE22) 


lie 


SIPPLE PIXEL. 


INPUT 

FVR 



R 2 

HI 

e.3 

FVR, THETA, WR 


99. 6e 

H2 

e.4 

FVR, THETA, WR 


99.21 

N3 

e.s 

FVR, THETA, WR 


99.46 

N4 

e.6 

FVR, THETA 


99.69 

H5 

e.7 

FVR, THETA 


99.79 

HI 

e.3 

FVR, THETA 


95. e4 

U 2 

e.4 

FVR , THETA , EFP 1 , WR 

- 

98. 6e 

U3 

e.s 

FVR, THETA 


95.19 

W4 

e.6 

FVR, THETA 


94.84 

U5 

e.7 

FVR, THETA 


97.98 

NS 

e.3 

FVR, WR,EFP1 


99. 7e 

N7 

e.4 

FVR, WR 


99.53 

NB 

e.s 

FVR, WR 


99.65 

N9 

e.6 

FVR 


99.67 

Hie 

e.7 

FVR 


99.75 

us 

e.3 

FVRjEFPl 


99.15 

W7 

e.4 

FVR 


98.81 

UB 

e.s 

FVR 


98.88 

W9 

e.6 

FVR 


99.47 

uie 

e.7 

FVR 


99.22 
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TABLE X- LONG. TOSILE STRENGTH (SCXOT) 


SIMPLE MODEL 

INPUT FVR itltfE AC mrlEl) 


R* 


N1 

0.3 

FVR 

N2 

0.4 

FVR, SFPT 

N3 

0.5 

FVR 

N4 

0.6 

FVR, SFPT, THETA 

N5 

0.7 

FVR 

W1 

0.3 

FVR, SFPT 

U2 

0.4 

SFPT, FVR 

W3 

0.5 

EFPljSFPT 

m 

0.6 

FVRjETP 

U5 

0.7 

FVR, SFPT 


N6 

0.3 

FVR, SFPT 

N7 

0.4 

FVR, SFPT, EFP1 

N8 

0.5 

FVR, SFPT 

N9 

0.6 

FVR, SFPT 

NIO 

0.7 

SFPT, FVR 

U6 

0.3 

SFPT, FVR 

U7 

0.4 

FVR 

U8 

O.S 

FVR, SFPT 

W9 

0.6 

FVR, SFPT 

M10 

0.7 

SFPT, FVR 


12.25 
43.72 
21.68 
43.68 
4®. 97 
33.37 
39.02 
26.13 
42.27 
33.55 

52.12 
68.43 
34.89 

49.00 

24.00 
46.61 
19.33 

33.13 
34.40 
37.65 
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TABLE XI 

- LONG. COCFRESSIVE STRENGTH (SCXHC) 


INPUT 

FVR 

SIMPLE MODEL 
TERMS ACCEPTED 

R* 

N1 

C.3 

FVR 

12.25 

N2 

0.4 

FVR 

18.23 

N3 

0.5 

NONE 


N4 

0.6 

SFPC 

8.52 

N5 

0.7 

FVR 

8.06 

W1 

0.3 

WR 

8.02 

W2 

0.4 

THETA 

9.29 

W3 

0.5 

GfF.STFC 

20.59 

m 

0.6 

THETA 

9.18 

U5 

0.7 

NONE 


N6 

0.3 

SFPC 

11.30 

N7 

0.4 

NONE 


N8 

0.5 

NONE 


N9 

0.6 

NONE 


NIC 

0.7 

GFP12 

12.01 

U6 

0.3 

FVR 

9.40 

V? 

0.4 

WR 

10.76 

UB 

0.5 

WR 

9.85 

W9 

0.6 

GFP12 

8.87 

W1C 

0.7 

NONE 
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TQPiir XII- TRANSVERSE TENSILE STRENGTH (SCYST) 


SirFLE ICDEL 



FVR 


R 2 

INPUT 

i£Rr£> 


N1 

0.3 

FVR 

27.03 

N2 

0.4 

FVR 

32.91 

N3 

0.5 

STPT 

8.10 

N4 

0.6 

FVR 

41.92 

NS 

0.7 

NONE 


Mi 

0.3 

FVRjWR, STPT 

26.89 

W2 

0.4 

FVR 

41.43 

U3 

0.5 

FVR 

14.74 

U4 

0.6 

NONE 


U5 

0.7 

FVR, bliP T 

31.05 

N6 

0.3 

FVR 

9.43 

N7 

0.4 

FVR 

8.19 

N8 

0.5 

FVR.EFP2 

15.58 

N9 

0.6 

NONE 


H10 

0.7 

NONE 


W6 

0.3 

FVR 

33.87 

U 7 

0.4 

FVR 

13.39 

MB 

0.5 

SKT 

8.62 

W9 

0.6 

FVR 

27.85 

U10 

0.7 

FVR 

32.77 
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TABLE XIII- TRANSVERSE COTFHESSIVE STRENGTH (SCYYC) 


SlfPLE KX)EL 


INPUT 

FVR 

TERTE ACCEPTED 

R a 

N1 

©.3 

FVRjSH’C 

33.17 

N2 

©.4 

FVR 

30.10 

N3 

0.5 

NONE 


N4 

©.6 

FVR 

38.93 

N5 

0.7 

NONE 


U1 

0.3 

FVR, WR 

28.19 

U2 

0.4 

FVR 

43.26 

U3 

0.5 

FVR.SrPC 

19.57 

U4 

0.6 

NONE 


U5 

0.7 

FVR 

15.85 

N6 

0.3 

NONE 


N7 

0.4 

NONE 


N8 

0.5 

NONE 


N9 

0.6 

NONE 


N10 

0.7 

HONE 


US 

0.3 

FVR 

28.68 

U7 

0.4 

FVR 

11.64 

US 

0.5 

NONE 


U9 

0.6 

FVR 

31.97 

Ul© 

0.7 

FVR 

33.05 
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TQHlg Xiv- in PLANE SHEAR STRENGTH (SCXY5) 


INPUT 

N1 

FUR 

0.3 

SirPLE MODEL 

teste accepted 

R 2 

FVR, THETA, OFP 12 

28.51 

N2 

0.4 

FUR 

8.74 

N3 

0.5 

THETA 

14.96 

N4 

0.6 

THETA, GFP12 , FUR , SIPS 

31.84 

N5 

0.7 

NONE 

• 

W1 

0.3 

THETA, UUR ,STPS, FUR 

48.16 

W2 

0.4 

FUR 

43.26 

U3 

0.5 

THETA 

8.40 

U4 

0.6 

THETA 

14.75 

U5 

0.7 

NONE 


N6 

0.3 

NONE 


N7 

0.4 

NONE 


N8 

0.5 

NONE 


N9 

0.6 

SIPS 

8.25 

N10 

0.7 

UUR 

8.53 

U6 

0.3 

FUR.STFS.OP 

29.06 

U7 

0.4 

NONE 


U8 

0.5 

NONE 


U9 

0.6 

GFP12.FUR 

22.20 

W10 

0.7 

STPS 

17.73 
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Further regression node Is were studied, in an attenpt to improve 
the predictive capability of the models, especially far the strengths. 
These models, incorporating higher order functions and combinations of 
predictor variables used in the simple models, show some improvement 
over the simple models, proving the value of including the "interaction" 
effects of predictor variables in the regression models. In addition, 
the higher order interaction models can fit response functions over a 
wider rarxje of fiber volume ratio, with associated improvements in the 
R* statistics. The data cases CONI and C0N2 contain selected points 
from the entire range of fiber wlune ratios, to supply the saitples for 
these runs. Furthermore, since higher order models are postulated, 

THETA is taken to be the cosine of the average of fiber orientation 
angles. The variable 1WR is a "dunny" variable, that is a function of 
other variables in the model. It is defined as 

WR = 1 - FVR - WR 

aid is intended to represent an "average" matrix voluie ratio over the 
thickness of the ply. The interaction models are shovn in Tables XV - 
XXVI. 

The general form of tte postulated models now includes higher order 
terms, so the predictor variables are tested up to the fourth pover. 
Symbolically, 

Y = B„ + B, (THETA) + B*(FVR) + B 3 (WR) ♦ Bq(EFPl) + B S (EIP) + 
B*(HVR) ♦ B 7 (THETA ) 2 ♦ B, (THETA) (FVR) + B, (THETA) (WR) ♦ 
B 10 (THETA) (EFP1) + ... + B,,(THEIA) 2 (FTO)(EFP1) + ... 

+ ... etc. 


BisflHETA)' 1 + B 13 (FVR) 4 
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’ll* number of terns possible in a complete fourth powr polynomial 
expansion becomes unwieldy for the cases studied. Considering the 
limitation of the size of the predictor matrix in the regression package 
used (100 x 100), the terms are intuitively grouped in the hope of 
eliminating large groups at one time. The regressions are conducted 
using “unlikely- candidates for admission into a particular model, and 
if no terms are entered, subsequent regressions are conducted without 
those terms. The justification for this approach is not a statistical 
argument, rather an interpretation of the physical principles active in 
any chosen model. The regressions to eliminate terms are merely used as 
a check on vhat seems intuitively reasonable. 
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TPlgg XV- LONGITUDINAL MODULUS (ECU) 


IKTERACTION MODEL 


INPUT 

FUR 

TERTE ACCEPTED 


R 2 

N1 

0.3 

THETA 14 *FVR*EFP 1 


84.50 

H2 

0.4 

THETA g *FVR*EFP 1 


92.66 

N3 

0.S 

THETA^FVRkEFP 1 


93.76 

N4 

0.6 

THETA a *FVR*EFP 1 


94.24 

NS 

0.7 

THETA‘ 4 *FVR*EFP 1 


85.08 

W1 

0.3 

THETA m *FVR*EFP 1 


63.84 

U2 

0.4 

THETA U *FVR*EFP 1 

* 

89.86 

U3 

0.5 

THETA' I *FVR»*EFP 1 


71.79 

m 

0.6 

THETA “ *FVR«EFP 1 


64.37 

U5 

0.7 

THETA “^FVRmEFP 1 


55.68 

N6 

0.3 

FVRhEFPI , EJf^HlWR 


99.82 

N7 

0.4 

FVR*EFP1 , FVR g 


99.83 

N8 

0.5 

FVR*EFP1 


99.72 

N9 

0.6 

FVRJtEFPl , EIf> 2 *VVR, WR g 


99.79 

N10 

0.7 

FVR«EFP1 , E1*»*WR 


99.79 

U6 

0.3 

fvr*efpi,wr 


99.17 

U 7 

0.4 

FVR*EFP1 , 1WR**FVR 


98.53 

UB 

0.5 

FVR*EFPi , JWR 2 *EH* 


98.99 

U9 

0.6 

FVR*EFP1 


99.58 

W10 

0.7 

FVRwEFPi , ETF*MVR 


99.38 

CONI 

VARIES 

THETA u *FVR*EFP 1 


96.48 

CON2 

VARIES 

FVR*EFP1,WR m 


99.92 
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Tpra f XVI- TRANSVERSE fCDULUS (EC22) 


INTERACTION IEDEL 


INPUT 

FVR 

TERfE ACCEPTED 

N1 

0.3 

FVR*EFP2*EIf > , EFP2 2 *FVR, THETA*FVR*IWR 

N2 

0.4 

FVR*EFP2*EIf , ,EFP2 2 *FVR, THETA 

N3 

0.5 

FVR*EFP2*ETF , THETA 2 *EFP2 , WR" 

N4 

0.6 

FVR*EFP2«ETP , THETA 2 *EFP2 , EFP2 2 *FVR 

N5 

0.7 

*** NEARLY SINGULAR 

W1 

0.3 

THETA 2 *MVR , EJF 2 *EFP2 , IWR 2 

U2 

0.4 

FVR*ETF , THETA , EPP2 2 *FVR 

W3 

0.5 

FVR*EFP2*ETF , THETA , FVR*EFP2«1 , WR 

W4 

0.6 

THETA 2 *WR , 1WR 2 *EFP2 , EIF* WR 

U5 

0.7 

*** NEARLY SINGULAR 

N6 

0.3 

FVR*EFP2*ETP , EFP2 2 *FVR , ETF*1WR 

N7 

0.4 

FVR*EFP2*ETF , ETP 2 FVR 

N8 

0.5 

FVR*EFP2*ETF , FVR*EFP2*IWR 

N9 

0.6 

FVR*EFP2*ETF , FVR^EJF 

NIO 

0.7 

FVR U , EFP2*ETF 

U6 

0.3 

FVR*EFP2*EIF , EFP2 , EFP2 2 «Elf J , FVR*EFP2 

W7 

0.4 

FVR#EFP2»ETF , FVR*EFP2*!WR , FVR«WR 

UB 

0.5 

FVR*EFP2wEJF , FVR#EFP2*WR 

U9 

0.6 

FVR«EFP2*E2'P t EFP 2 a *rKR 

M10 

0.7 

FVR*EFP2*EMP , WR 2 *FVR 

CONI 

VARIES 

m NEARLY SINGULAR 

COM2 

VARIES 

FVR*EFP2*EJP , FVR*EFP2*IWR 


R 2 


99.19 

99.55 

98.92 

99.22 

93.26 

96.79 

93.49 

88.35 


99.22 
99.07 
98.89 
99.14 

99.23 
98.62 
98.28 
97.93 
98.44 
97.86 


99.79 
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TABLE XVII- IN PLANE SHEAR MODULUS (EC12) 


INTERACTION MODEL. 


INPUT 

FVR 

HiRfE ACCfcKlED 

n 

N1 

0.3 

THETA , FVRKaT , THETA 1, *FVR*<3f > 

97.86 

N2 

0.4 

*** NEARLY SINGULAR 


N3 

0.5 

THETA 1 , FVR 2 *GW 5 f CFPi2*GrF 

97.75 

N4 

0.6 

THETA, FVR 2 *WR 

98.01 

N5 

0.7 

THETA, FVR 2 «<3 , F,GFP 12 

98.46 

HI 

0.3 

THETA' 1 , FVR‘ 4 *K3TP,FVR 2 

95.49 

W2 

0.4 

*** NEARLY SINGULAR 


U3 

0.5 

THETA U *FVR*IWR,THETA U *WR, WRKGfP 

91.04 

m 

0.6 

THETA u , FVR 3 *<3T , THETA 

96.70 

U5 

0.7 

*** NEARLY SINGULAR 


N6 

0.3 

FVR*GMP , WR 2 «GFP 1 2 

97.73 

N7 

0.4 

FVR*GfF,GFP12 

97.97 

N8 

0.5 

FVR*GtF , FVR*KTP12 

96.52 

N9 

0.6 

FVR*Q1P,FVR 2 *<3 : Pi2 

97.10 

N10 

0.7 

FVR’ , «Grt > ,FVR«GFP12«GW,WR*«GW 

98.90 

US 

0.3 

FVR«GMP , FVR**GFP12 

96.91 

W7 

0.4 

FVR*GTP , FVR*GFP12 

92.37 

UB 

0.5 

FVR^GfP , FVR*GFP12 

95.08 

W9 

0.6 

FVR#GflP,FVR 2 «GFP12 

97.42 

HI® 

0.7 

FVR a *K»,FVR»*(FP12**<M > ,IWR 2 «<a*» 

96.85 

CONI 

VARIES 

FVR 2 *WR, WR*WVR , FVR#GTP , THETA' , *FVR#«f J 

99.09 

C0N2 

VARIES 

FVR 2 «OT , WR*a*> , GFP1 2*K3f > 

99.54 
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TPFgF XVIII- LONG. THEKHQL BffflHSIOH fCTCll) 
INTERACTION 1CCEL 


INPUT 

FVH 

TERTB ACCEPTED 

R 2 


N1 

0.3 

THETA 2 *MVR , WR 2 , FVR*ETP1*!WR, E2f» 2 *EPPl 

92.51 

K2 

0.4 

THETA 2 ^WR,THETA g l EFPl 1, ,Erff ,2 «FVR,fWK 2 »EFPl 96.38 

N3 

O.S 

THETA 2 *rWR f WR,Q* ,a,( WH|E* , * ,2,tE3:pA 

97.26 

N4 

0.6 

FVR*EFP1 , THETA*FVR»EFP 1 ,EMP 2 *WR 

96.32 

H5 

0.7 

THETA, ETF*MVR 

90.66 

W1 

0.3 

*** NEARLY SINGULAR 


U2 

0.4 

THETA", rWR 2 «EIF 

80.81 

U3 

0.5 

theta", ejf 2 *iwr, 

87.98 

m 

0.6 

theta" ,m 2 *wr 

75.20 

W5 

0.7 

THETA", FVR 2 *IWR 

82.97 

N6 

0.3 

MyR 2 «ETP , EPP1 ^EM* , FVR" 

99.29 

H7 

0.4 

HVR 2 «ETP , FVR*EFP1*MVR , IWR 2 *WR 

99.17 

N8 

0.5 

wr 2 *etp,eppi,fvr" 

98.94 

N9 

0.6 

mur 2 *etp , EFP1 2 *rWR 

98.94 

N10 

0.7 

mur 2 «emp , EFPi*rwR • 

99.33 

U6 

0.3 

MVR 2 *EIP,FVR*EFPl*tNR 

98.35 

W7 

0.4 

HVR 2 «ErF , FVR*EFP1«MVR 

98.55 

US 

0.5 

HyR 2 *ETP , EFP1 2 «JWR 

98.56 

W9 

0.6 

MVR 2 *EJF,EFP1 2 *1WR 

99.00 

W10 

0.7 

MVR 2 *EFP1 ,EMP 2 *MVR 

98.20 

CONI 

VARIES 

THETA, HUH 3 ,EFP1 2 *WR,FVR*WR*EFP1 

96.82 

CON2 

VARIES 

MVR 2 *EIP,FVR*EFP1»IWR,FVR" ,FVR 2 *WR. . . 

99.84 
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TOCTF vtx- TEAMS. TWEHHQL EXP AMS I OH (CTE22) 


INTERACTION IC C EL 


INPUT 

FVR 

TERTE ACCEPTED 

n 

N1 

0.3 

theta 2 *!*®, wr 

99.60 

N2 

0.4 

THET , A I *MUR,WR 2 *FVR J FVR»EFPl*EIf > 

99.38 

N3 

0.5 

theta 2 *!*®, mur 2 *fvr, fve" 

99.48 

N4 

0.6 

THETA 2 *!*® , !*® 2 *FVR , EMP **VUR, I*® 2 *WR 

99.73 

N5 

0.7 

FVR 2 , THETA , THETA 2 *FVR 

99.81 

W1 

0.3 

theta 2 *nyR , UWR , ETP 1 2 *wr 

95.16 

W2 

0.4 

theta 2 *!*® , I*® , FVR*EFP1*!*® 

98.71 

U3 

0.5 

THETA 2 *!*® , I*® 2 *FVR, EPP 1 **WR 

95.91 

W4 

0.6 

THETA 2 *!*®, THETA, THETA* 1 

95.69 

W5 

0.7 

*** NEARLY SINGULAR 


N6 

0.3 

FVR,!*® 3 

99.70 

K7 

0.4 

FVR, I*®* 4 

99.59 

N8 

0.5 

FVR,!*® 2 

99.67 

N9 

0.6 

FVR,EMP 2 *EFP1 

99.70 

N10 

0.7 

FVR, l*®* 4 

99.82 

U6 

0.3 

FVR, FVR*EFP 1*EW 3 

99.26 

W7 

0.4 

FVR,!*®* 4 

98.97 

UB 

0.5 

FVR 

98.88 

119 

0.6 

FVR,!*® 14 

99.57 

W10 

0.7 

fvr,etp 2 *fvr 

99.29 

CONI 

VARIES 

THETA 2 *!*® 

99.32 

C0N2 

VARIES 

fvr,fvr 3 ,i*® 2 *etp 

99.95 
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TABLE XX- POISSON RATIO: MAJOR FNUC12) 


INTERACTION MODEL 


INPUT 

FVR 

TERTB ACCEPTED 

R* 

N1 

0.3 

*** NEARLY SINGULAR 


N2 

0.4 

THETA , EFP2*MVR 

97.96 

N3 

0.5 

THETA, GPP1 12*1WR 

96.71 

N4 

0.6 

THETA, EFPl*l¥JR 

98.17 

N5 

0.7 

. THETA, FVR*EFP2 

96.48 

W1 

0.3 

*** NEARLY SINGULAR 


U2 

0.4 

THETA , THETA' 4 *FVR*GFP12 

84.73 

U3 

0.5 

THETA 

89.43 

m 

0.6 

THETA , WR*GFP 1 2 

84.27 

U5 

0.7 

THETA , FVR*MWR 

92.10 

N6 

0.3 

FVR 

97.83 

N7 

0.4 

FVR 

98.48 

N8 

0.5 

FVR 

97.77 

N9 

0.6 

FVR,FVR#1WR 

98.52 

NIO 

0.7 

FVR 

99.17 

MS 

0.3 

FVR 

97.32 

U 7 

0.4 

FVR,WR»EFP2 

96.50 

UB 

0.5 

FVR 

96.38 

W9 

0.6 

FVR , EFP i*ETP2 , GFP 1 2*IWR 

98.41 

U10 

0.7 

FVR,FVR*EFP2 

96.97 

CONI 

VARIES 

mw SINGULAR 


C0N2 

VARIES 

fWR, FVR*1WR,EFP1*MVR 

99.77 


TABLE XXI- POISSON RATIO; MINOR (NUC21) 


INPUT 

N1 

FVR 

0.3 

INTERACTION MODEL 
terte ACCEPTED 

R* 

THETA , THETA “"FVRkEFP 1 

91.69 

K2 

0.4 

THETA , FVR*EFP 1 

94.66 

N3 

0.5 

THETA, FVR*ETP1 ,EFP2*GFP12 

95.10 

N4 

0.6 

THETA , THETA a *FVR*EFP 1 ,EFP2 

97.15 

N5 

0.7 

THETA, THETA a *FVR*EFP 1 

95.82 

W1 

0.3 

THETA, FVR*<3FP 12 

91.16 

M2 

0.4 

THETA , EFP2*WR 

89.52 

M3 

0.5 

THETA , FVR*EFP1 

92.06 

M4 

0.6 

THETA 

92.53 

M5 

0.7 

THETA , FVR*EFP 1 , THETA 14 *FVR*WR , EFP 2* WR 

95.60 

N6 

0.3 

FVR*EFP1 , FVR*EFP2 

95.48 

N7 

0.4 

FVR*EFF1 , FVR*EFP2 

94.69 

N8 

0.5 

FVR*EFP1 , FVR*EFP2 , FVR*GFP12 

95.52 

N9 

0.6 

FVR*EFP1 , FVR*EFP2 

92.85 

N10 

0.7 

FVR* EFP 1 , FVR*EFP2 

91.77 

W6 

0.3 

FVR*EFP1 ,GFP12*MVR 

87.83 

M7 

0.4 

FVR*EFP 1 , FVR*EFP2 , FVR*rWR 

86.48 

MB 

0.5 

FVR*EFP1,EFP2 

84.36 

M9 

0.6 

FVR*EFP1 , FVR*EFP2 

89.84 

M10 

0.7 

FVR*EFP1 , FVR*EFP2 

91.55 

CONI 

VARIES 

THETA,FHV*GFP12,EFF2 f THETA ,, *FUR*IKR, . - . 

98.70 

COM2 

VARIES 

FVR*EFPl,FVR*rWR,EFP2,WR*GFP12 

98.35 
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TABLE XXII 

- LONGITUDINAL TENSILE STRENGTH (SCXXT) 




INTERACTION MODEL 


INPUT 

FVR 

TERMS ACCEPTED 

R 2 

N1 

0.3 

THETA a *FVR»SFPT 

17.72 

N2 

0.4 

THETA 14 *FVK*SFPT , MUR* 4 

47.65 

N3 

0.5 

THETA a *FVR*SFPT 

27.65 

N4 

0.6 

THETA l4 *FVR*SFPT , FVR 2 *EMP 

44.67 

N5 

0.7 

THETA M *FVR*SFPT 

45.35 

W1 

0.3 

FVR*SFFT , FVR#EFP 1*1VR 

39.18 

U2 

0.4 

FVR*SFPT , FVR 2 *IWR 

42.87 

W3 

0.5 

EFP 1*SFPT, EJF«IWR 

33.97 

m 

0.6 

EIF 2 *IWR , THETA 2 *SFPT»WR 

45.09 

U5 

0.7 

theta 2 *fvr*sfpt 

32.56 

N6 

0.3 

FVR*SFPT , FVR 2 *EJP 

52.95 

N7 

0.4 

FVR«SFPT , FVR*EFP 1 

64.41 

N8 

0.5 

FVR*SFPT, WR 2 *FVR 

39.12 

N9 

0.6 

FVR*SFPT 

47.13 

N10 

0.7 

FVR*SFPT, FVR*EMP 

27.43 

W6 

0.3 

FVR«SFPT , WR 2 *EFP 1 

49.71 

W7 

0.4 

FVR*SFPT 

25.19 

U8 

0.5 

FVR*SFPT 

32.16 

U9 

0.6 

FVR*SFPT 

34.06 

W10 

0.7 

FVR*SFPT 

35.09 

CONI 

VARIES 

THETA w *FW*SFPT f FVR*VVR*1WR, FVR*EFP1*WR 

81.20 

COM2 

VARIES 

FVR*SFPT, FVR»WR f WR 2 *SFPT 

84.79 
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TQPiir XXIII- LONGITUDINAL OUtfRESSIVE STRENGTH (SCXXC) 


INTERACTION MDDEL 


INPUT 

FVR 

TERIE ACCEPTED 


R a 

N1 

0.3 

SFPC*SWC 


12.53 

N2 

0.4 

FVR«1WR 


19.45 

N3 

0.5 

NONE 


9.81 

N4 

0.6 

SFPOQT 


N5 

0.7 

FVR<KTP12 


10.20 

U1 

0.3 

wr" 


10.40 

U2 

0.4 

theta 14 

- 

9.32 

U3 

0.5 

ETFkSITC , THETA**EMP 


23.32 

U4 

0.6 

THETA' 4 


9.20 

U5 

0.7 

NONE 




N6 

0.3 

SFPC , C3T>12*SrPC 

20.04 

N7 

0.4 

NONE 
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0.5 

NONE 


N9 

0.6 

NONE 


N10 

0.7 

GFP12*EMP 

14.96 

U6 

0.3 

FVR*MVR 

11.91 

117 

0.4 

WR 

10.76 

UB 

0.5 

WR 

9.85 

U9 

0.6 

CTP12 a 

9.10 

H10 

0.7 

NONE 


CONI 

VARIES 

FVR*WR 

46.48 

C0N2 

VARIES 

FVR* WR , SFPC 

44.44 
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tout XXIV- TRANSVERSE TOSILE STRENGTH (SCYYTl 




INTERACTION IEOBL 


INPUT 

FUR 

TERTE ACCEPTED 

R 2 


N1 

0.3 

WR 2 «SPFT 

31.60 

N2 

0.4 

HVR 2 *SfFT 

37.23 

N3 

0.5 

em^sift 

9.61 

N4 

0.6 

fvr 2 *iwr 

47.59 

N5 

0.7 

NONE 


U1 

0.3 

FVR 2 «WR,SW r r 

25.39 

W2 

0.4 

SMPT*WR 

43.94 

U3 

0.5 

wr 2 *fvr 

16.32 

U4 

0.6 

FVR*WR*EFP2 , ETP 

24.10 

U5 

0.7 

FVRMSTPT 

30.29 

N6 

0.3 

FVR«rWR 

10.47 

N7 

0.4 

wr 2 *fvr 

8.94 

N8 

0.5 

FVRwEFP2*SITT 

13.54 

N9 

0.6 

sift 2 *muk 

9.40 

N10 

0.7 

FVR**ETP 

9.13 

U6 

0.3 

jfr 2 *eip 

35.13 

W7 

0.4 

mur 2 *wr 

19.34 

UB 

0.5 

SIFT 2 *!*® 

12.89 

W9 

0.6 

FVR 2 «EFP2 

29.27 

U10 

0.7 

wr 2 «sift 

36.77 

CONI 

VARIES 

theta 11 «sm>WR , fvr*efp2*wr, *srFT J «rwR 

73.42 

CON2 

VARIES 

SfFT 2 *lWR , FVR*WR*IWR 

76.40 
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TflaF xxv- TRflhsvfeflSE COryiiEsSIVE siHE NulH (SCYYC) 
INTERACTION YODEL 


INPUT 

FVR 

tertb ACCEPTED 

XV 

N1 

0.3 

SfPC*IWR 

33.39 

N2 

0.4 

fvr 2 *etf 

32.99 

N3 

0.5 

NONE 


N4 

0.6 

FVR**IWR 

42.31 

N5 

0.7 

NONE 


U1 

0.3 

FVR 2 *WR 

26.24 

U2 

0.4 

fvr^ejf 

43.86 

U3 

0.5 

S«PC**WR 

21.13 

m 

0.6 

FVR*WR*EFP2 , ETP 

25.75 

U5 

0.7 

SMPC*IWR 

18.63 


N6 

0.3 

Slf»C a «WR 

11.57 

N7 

0.4 

EFP2*IWR 

9.03 

N8 

0.5 

FVR*EFP2 

9.87 

N9 

0.6 

NONE 


N10 

0.7 

SfPC^lWR, FVR a *IWR 

19.07 

U6 

0.3 

wr 2 *etp 

32.50 

W7 

0.4 

EFP2*IWR 

14.58 

UB 

0.5 

NONE 


U9 

0.6 

fWR 1 «SK»C 

32.85 




35.79 

U10 

0.7 

wr 


can 

VARIES 

THETTA U *SI , FC*WR , FVR* 1 

76.43 

C0K2 

VARIES 

MVR J , FVR a *lWR, IWR^SWC 

75.59 
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TABLE XXVI- IM PLANE SHEAR STENGTH (SCXVS) 


INTERACTION MODEL 


INPUT 

FVR 

TERfS ACCtyiED 

n 

N1 

0.3 

FVRwGFPl 2#GTF , THETA 1 * 

27.64 

N2 

0.4 

FVR#CFP12*EMP 

13.51 

N3 

0.5 

THETA 

14.97 

N4 

0.6 

THETA l| «<3 r P12,Sri , S*rWR 

30.84 

N5 

0.7 

NONE 


U1 

0.3 

THETA, FVR# WR#EMP, THETA l, #SH , S,FVR#IWR 

52.20 

U2 

0.4 

THETA 14 #FVR , THETA g #OFP 1 2 

26.58 

W3 

0.5 

THETA* 4 

12.89 

W4 

0.6 

THETA, FVR# WR 

22.33 

U5 

0.7 

THETA 14 #FVR 

10.72 

N6 

0.3 

NONE 


N7 

0.4 

Sff*S#WR 

11.24 

N8 

0.5 

NONE 


N9 

0.6 

si-fs, sips'* 

16.14 

NIO 

0.7 

FVR 2 «WR 

11.40 

U6 

0.3 

srps#wR l Grp‘ l 

28.58 

U7 

0.4 

FVR#<FP12#MVR 

8.28 

U8 

0.5 

NONE 


U9 

0.6 

FVR«GFP12*MWR 

19.20 

W10 

0.7 

sjps 

17.73 

CONI 

VARIES 

THETA * , #FVR , FVR 2 #SIPS 

36.74 

C0N2 

VARIES 

FVR«WR, lWR a , FVR 2 »WR 

61.46 



CHAPTER IV 


DISCUSSION 


A. Overview 

The miner i cal simulations conducted show that certain assumptions 
about the statistical distribution of local nonuniformities in fiber 
conposites lead directly to quantifiable variations in material 
properties. The advantages inherent in the stochastic characterization 
are nunerous. The development of quality control and reliabilty 
measures for conposites is crucial to their acceptance in aircraft 
designs. The reduction in needed experimental data achievable through 
judicious simulation of the wide variety of available composite material 
systems could significantly lower the costs of material selection and 
acceptarce testing. In the results of this study, the confidence 
intervals calculated can be interpreted as the product of an 
experimental program, specifically designed as an analog of the physical 
processes which occur in real materials. 
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B. Histograms and Distributions 

Data cases 1, 2, and 3 demonstrate the differences between a 
deterministic base case and random cases with narrow and wide dispersion 

of input data about the base case. 

In Fig. 30, it is apparent that the deterministic case 1 wall® of 
1575© ksi. for longitudinal nodulus falls near the mean of the case 2 
data. However, the case 3 sample appears to hax® a nean slightly lover 
(approximately 15000 ksi.). It should be noted that the size of the 
interval over v* ich the sample occurs is noticeably larger in the widely 
distributed case 3 run. 

Transverse nodulus, (Fig. 31) demonstrates a higher mean valt® for 
the wide distribution than for the narrow, vhich is greater than the 
deterministic valt® of 1065 ksi. reported in Table II. The increased 
transverse modulus is related to the added stiffness available in fibers 
with high misalignment relative to longitudinal direction. 

Shear modulus, (Fig. 32) is measurably changed by nonuniformities. 
The deterministic value of 516 ksi is exceeded by the case 2 value of 
approximately 620 ksi, viiich is further exceeded by the case 3 vail® 
near 900 ksi. Fiber misaligment has a significant effect in shear 
nodulus variation. 

Poisson's ratios (Fig. 33, 34) show similar trends in location of 
sample means and relative dispersion of the sample for the data studied. 
Poisson's ratios generally increase with fiber misaligment and volume 


fraction changes. 
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The coefficients of thermal expansion (Figs. 35, 36) for the sample 
studied reflect the longitudinal contraction of graphite fibers tdien 
heated. The longitudinal coefficient of thermal expansion for 
AS-graphite fiber is -0.550 x 10 _6 / F, while the transverse coefficient 
is 0.560 x 10~ 5 / F. The offset orientation of crystal lattice planes 
in graphite fibers can explain this behavior. These values, the fiber 
misalignment, and fiber volume ratio near 0.5 all contribute to the 
occurrence of a negative longitudinal coefficient of thermal expansion 
for the composite, fit higher fiber volume ratios, the values calculated 
would be less than in the present case, because of the contolling fiber 
behavior for high fiber volume ratio. 

The longitudinal strengths (Fig. 38, 39) are significantly reduced 
when nonuniformities are present. The deterministic case 1 value of 203 
ksi. for tensile strength is compared to a mean near 160 ksi for case 2 
and a nean near 130 for case 3. In compression, the deterministic value 
of 165 ksi. compares to means near 100 ksi. and 80 ksi. for the narrow 
and wide distributions, respectively. The failure mode in compression 
varies in the random sanples. 

Transverse strengths (Fig. 40, 41) show sensitivity to the 
variations assumed. Misalignments, volume fraction nommiformities, and 
constituent strength variations all contibute to reduction in the 
strength values. Sid) -ply shear failures occur, which undermine the 
already low transwrse composite strengths. 

In plane shear strength (Fig. 42) values decline from 10.01 ksi. 
for case 1 to a mean near 8.0 ksi. for case 2. However, case 3 shows a 
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value of a nean near 8.0 also. It appears that the added shear strength 
due to fiber misalignnent is balanced by the reduced strength d«s to 
variable fiber wlune fraction. 

C. Confidence Curves 

effects of various shape parameters of fiber strength are shovm 
in Figs. 43 and 44. Tte higher weibull distribution shape parameter of 
20 produces a narrow distribution of fiber strength values. The 
conposite that has few weaker fibers is expected to he stronger, and 
Fig. 43 demonstrates this for lonitudinal tensile strength. Hovever, 
conpressive failure (Fig. 44) is a nore complex phenomenon. In the 
region of low fiber volume ratio, the ’rule of mixtures’ failure 
criteria for a subply can control the failure mode. At higher fiber 
volune ratio, however, conpressive failure can be initiated by 
delaminat ion , or by a shear failure in a sub-ply. The mixture of 
failure modes in conpressive failure is not well understood, but can 
explain the seeming inconsistency of the intersection of the curves in 
Fig. 44. At a fiber volune of 0.7, the weakest fibers (a = 10) are in 
the strongest conposite, %*en strength is normalized with respect to 
fiber conpressive strength. 

Tte effects of various shape paraneters for iratrix strengths are 
studied in Figs. 45, 46, and 47. Transverse tensile and conpressive 
strengths show expected reductions for lowr matrix strengths. In-plane 
stear strength shows lower dispersion at a large fiber rolune of 0.7, 
and also declirms in general for higher fiber volune. 
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The fiber misalignnent effects are stulied in Figs. 48-57. 
Longitudinal modulus (Fig. 48) shows narrow intervals and slight 
reductions for greater misalignnent. Transverse modulus (Fig. 49) and 
in plane shear modulus (Fig. 50) are enhanced by fiber misalignnent. 
Longitudinal tensile arm! compressive strengths are degraded by 
misalignment (Figs. 51, 52). Transverse tensile and compressive 
strengths are enhanced (Figs. 53, 54). In-plane shear strength shows 
total separation of confidence intervals betveen curves with different 
degrees of misalignment. Poisson’s ratios (Figs. 56, 57) increase for 
high fiber misalignnent values. 

The fiber stiffness effects (Figs. 58—67) are wry smell for the 
distribution parameters studied. 

D. Examination of Regression ffadels 

The regression models for thermoelastic properties demonstrate 
resonably high predictive capability in the simple models assumed. 
Marginal improvements are achieved in expanding the models to include 
higher order interaction terms. Further improvement is gained by using 
sample data from the wide range of volume percent valiss. The higher 
multiple correlation coefficients of these models nay be due to the 
additional information available in the sanple size of 100 that mas 
uset j m The nearly singular predictor matrices thich occur in the higher 
order models indicate that terms must by selectively removed to 
eliminate linearity between assumed predictor terms. The regression 
results support the use of the simple models for themoelastic 
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properties, because improvements in predictive capability in the higher 
order node la for the same data are small. 

Strengths are not node led well by the simple or the interaction 
iiodels. The predictors chosen are average properties, vhereas the 
strengths are based on the weakest points in the material. Even the 
unidirectional cases (H6-N10, W6-W10) present data that the interaction 
nodels have considerable difficulty in accomodating. Somewhat greater 
predictive value is gaimd by using the expanded data for strength model 
prediction. Using fourth order algebraic functions, values of the 
multiple correlation coefficient square approach 857. for longitudinal 
tensile strength. The otter strengths generally have poorer results. 



CHAPTER V 


CONCLUSIONS 


A tractable, constituent based, probabilistic analysis procedure 
for fiber composites has been developed using the ICAN program as a 
basis. Uithin the limitations of the mechanics of material model, 
properties and strengths of a variety of composite material 
configurations can be simulated. 

This study quantifies the thernoelast ic and strength properties of 
a graphite/epoxy ply subject to assumed uncertainties for fiber 
misalignment, constituent volume fractions, and constituent properties. 
The results show several advantages of probabilistic characterization of 
this material. These include the identification of unforseen variations 
in composite material properties, and the mechanical effects of local 
nonuniformities. The relative importance of the various fabrication and 
material variables on composite properties is identified, and the 
resulting behavior quantified. 

The advantages of a probabilistic formulation of composite material 
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properties over a deterministic one are nunerous. Coiiparison of the 
results of this study with test data oould reveal sone souroes of 
previously unaccounted scatter in the data. Expected value ranges could 
be predicted for experinental results. Since the simulations provide 
data that is analagous to experinental data at lover cost, laboratory 
classification, material selection, and acceptance testing of composites 
can be guided by tte information made available by these methods. 

Although the method presented provides results for only the basic 
ply, extension of the simulation to incline lamination angle variations 
in a general layup is feasible. Since finite element material property 
cards are generated, structural analysis of components with randomly 
varied properties defined at a number of points in the body can supply a 
mure realistic description of the random nature of striutural response 
d«? to material inhomogeneity. 

Tfce stochastic formulation of material properties is generally 
recognized as ow» requirement of failure theories for materials. 

Although tte failure criteria in the models used in this stixiy are 
conservative, progressive failure of fiber composites could be modeled 
by incorporating load redistribution and material property recalculation 
in the vicinity of failed material. 
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nooooonnonoononoonn 


„ - 

"rriTiTir? 

PROBABILISTIC INTEGRATED COMPOSITES ANAL«ER <_P I 

COMP^ESSM^ 

?^S2| TE THE ANALYSIS SAMPLES TROM INPUT DISTRIBUTIONS TO OBTAIN 
COHPOSiTE^PROPERTIES AND GEOMETRY < WHICH ARE THEN INPUT TO ICAN. 

»u?r^ni|TpUT^INCLUDES^OUTPUT^DATASETS Or ICAN WHICH ARE NAMED 


M M 
MU 
M M 
M M 
MM 
MM 
MM 
MM 

M M N M M M M M M N N M M 
kmmmmmmmmmmmmmmmmmmmmm 


VMT9 M ft MASTER PROGRAM rot "ICAN" WHICH ALLOCATES 
DYNAMICALLY SUFFICIENT STORAGE FOR THE ARRAY VARIABLES 
IH "ICAN" AND "PIC AH" CODES. 

SS55 

I1AXLEH • MOO 
CALL SPINIT 
STOP 



C mD S IKPUT T DmiET M tO DETERMINE tr RROIAIILISTIC ANALYSIS IS DESIRED 
common /uin«/ inntdi . ourr . iwr . «r r . inds . idbk 

INTE0E« C *IH . RUNS # RUNMO.OUTf 

S \lin.Yuu cwm 

READ (SIN. USD .Ht.HUnS.SSMT 

ir (.NOT. BSTAT) 00 TO SSS 

C SERB (HH.ISUI RUNS 
C SET UP POINTERS roi RASTER RRRRT 

till ■ UI) • NUNS 
LISI • t(2> ♦ NW« 

UR) •1(1) ♦ NUNS 

US) > UR) • HUMS 
LUI • US) ♦ NUNS 
lin • UR) * HL 
us) ■ un * Nt 
u • u» 

tt • ui) 

LI • US) 

L« • UR) 
ts • US) 

M - UR) 

t? • UI) 

C LOOS* 'SUNS' TINES THROWN DATA CSERTION AND ICAN ROUTINE 

"KttWm!i?lR.MW».SIM).MtS).MMI.S(in.S*UI.Nl.NUW) 
REWIND 10BK 
CALL ICANNN 
ENDriU OUTf 
III CONTINUE 
REWIND ourr 
CO TO ARSt 
SIS CONTINUE 

CALL ICANNN 

:::: kkk 

|M) rORNAT(AR) 

ARM CONTINUE 
RETURN 
END 
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.ssss 

1NTE0ER ML.MLC.MrtS.lHT.l*. INPf 
1NTE0ER PIN . ““{JlS 110 

ss&wigr 


c 

c 


if ^B in unworn rrndoh nuhrer gwemioi 
remind iseedp 

RERDt ISEEDr.*) vm 
REMIND INPf 

M-jraSiKIiS!,.. 

RCRDCPIN.V) IDENTUl.Nt.NlCfNrtS 
irlNL.ED.NH5) CO TO J* 

MRZTC ClNPf.M* 

STOP 

II MRXTE CINPf.lt) IDENTCII.Nt.NLC.NHS 

RERDCPIN. ID COHSRT.MQLW.fMWV.VORRTV.COHV 

lUJJJriNUlI^CSRHIuTNHM.TICSIO.VfPHU.VrPSIQ. VVPLRH.KVVP 
MRITt CINPf.ll) CSRNR 

HEAD IflHill) MJJ 

uhitc (iHpr.m hw 
KCID (MHiin IIHDV 
MR1TCC INPf . I DJggy 

»T&»!.WS«f 

*RERD^CPlH.lD*MEMTCX).TU.TCU.PtC)l» D.®®* * ' 

If IIDENTCI) .M.PLT) CO TO R» 

WRITE CINPf »t! IDENTCJ) 

STOP 

! M ir IRHOIEV) 00 TO RR 
DO IMR • l.HL 
THETRCI*) • THnU 
•X CONTINUE 


CM 



CO TO HI 

ga !ssji”».m>».i«K-*> 

theta < i> > • » 

>•> ) THETM nUMl.l 

III CONTINUE 

c ir iriMTUi oo to ii* 

DO III II ■ l*"” 5 .. 
vrm«i • vrpnu 

Ul CONTINUE 

60 TO 1W ^ 

11% DO U9 !• ■ t*Nt1S 

115 cut UAANDCXll VPPSIO Yl 

ir u xt. •.»•» c® T0 1,5 
.vrrua) • t 

129 CONTINUE 

c no ir (vodatyico TO !*• 

DO 119 II ■ !••** 
warn) • VVPUH 

»• gran. 

' “* Sll"»WrU!5w;.«f<l«>> 

VVP(IA) • WPCIAI/HA* 

!«• CONTINUE 

' “,a8®i!fcABaaw«»fi‘ 

2*1 CONTINUE 

C Kll LOADING CONDITIONS 

°°»EH IPW. !i| t IOENTI JI.NDSn .«> I* »*? m!s* 
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no non 


JH CONTINUE 

hem output options 


lUimNitl) IDENTIC). IOUT 
UNITE <IHRr.II> IDENT(S)»IOUT 

INCREMENT MID REriLE SEED TOR TUTURE RUNS 

ISEED * ISCED ♦ 111 

REMIND ISEEDF 

URITE (ISEEDr.l) ISEEO 

2 FORMAT (IX.IIHIDENT(2) •.*•> 
i FORMAT <Ii> . 

• format t •'there is » nix ur in the layer properties card- > 

i ropnAT tAi.iii) 

II FORMAT (AD. Ill > 
it roinAT <u> 

!5 BBS? !u:S:»}S.n i i.j<i».n.».a.n- 

!> BBS? B. 

u rORMAT <A».2At, UX.2AA. jri.2> 

II roinAT <AI.2AI.2ri.2.2AA« sri.2) 

ii roinAT ui.iri.i) 

II FORMAT 
21 rORMAT (AI.IIJ 

21 rORHAT (AI.II) 

22 FORMAT (III 

\i roinAT "‘input ERROR. . . nms must be set equal to nl. • > 

RETURN 
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SUfeXOUTlME UMND(2{ hommmnmoniimiimiimmmhn 

om i/i/ 

INTEGER ft. X ... 

COMMON /SEE®/ iseed 

ir « .n. •> oo to i 

X • • 

n ■ i"M 

rn • n 

X • ISEED 

ft ■ ♦ 3 


I « - noDU'x, m 
rx • x 
z m rx/m 

RETURN 

END 


9&T 



n ftnnfl 



urn ri.nu.sictu.xi.xi.v 

yV»si«wM i'2«MOO<xi i »«• »>•<««< 2 -pi-xi ) ) ) *nu 

RETURN 

END 
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SUBROUTINE ORN(RtRNDR.K.K) 

C SUBROUTINE TOR OENERRTINO ORMMR BRRIRTES WITH PRRRMETERS 
C BLRtIDR RND K. 

DIMENSION Ut lit) 

DIMENSION HIM) 

COMMON /SEED/ I3ECB 
DO »B I • ).K _ 

91 CRLL URRHD(Ud) ) 

Ml) ■ U(l> 

DO ill I • 2.K 

IBB Ml) » U(I) ■ PII-I) 

* * • (-1 .B/RIRMDR) * RtOOCMK)) 

RETURN 
END 
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SUBROUTINE WEIBt XI . ALPHA. BETA . Y I 


_ 

i 

C VARIABLE DISCAIPTIOHS 
C ALPHA » SHAPE PARAMETER 

c ,C *| '• UNirOBntt*DISTBIBUTEO RANDOM VARIABLE ON (•■>) 

C I • HEIBULL DISTRIBUTED RANDOM VARIABLE 

C USE IS MADE OT THE HEIBULL DISTRIBUTION TUNCTION 

C F«) ■ 1 - EXPf - (X/BETA) »" ALPHA I EOR X .OE. EERO 

c 

t^BETA i I'-AUHItOMXII) -» H/ALPHAI 

RETURN 

END 
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SUBROUTINE COM 

»* .««» »«* ™* r,it ,o “ 


THIS ROUTINE SinPtT COPIES TUfc inr«» 


dimension wentisi 

too ICB L E «*NB?COnS»T^BIDE. BIHDW.HOMUDr 
INTEGER HL.NLC.NftS. INT.IR 
INTEGER BIN, POUT 
Dm PIN/SI/. POUT/5/ 

REftt TU.TCU.NBS.I«S 

Dm PET/* H»v 

READ! PIN. A) (DECK(II«I*1 ,20 . 

WRITE (POUT. I) (#KK(IM , l»n) 

BEBD(PiN.t) nns 

WRITE (POUT. 10) IDENT( 1 I .HL.MLC.nrn 

REAO(PIH.ll) C OHS IT 
WRITE (POUT. >J> COWSRT 
BEAD (PIN. It) CSBNB 
WRITE (POUT.ISI CSRNB 
READ (PIM.lt)BIDE 
WRITE (POUT. IS ) RIP E 
READ I PIN. IS) RINDV 
WRITE! POUT • I 1 IRINOV 
READ! PIN. 12 I NONUOr 
WRITE (POUT, IS) NOHUDf 

READ UTER DATA 

IBB READ (PIN, tO I0ENT(2l.INP(IR).IDnR).TU.TCU.PtO2.IR).tNETA(IBI. 

‘l“ (MENTIS). NE.PMl UD » IBS 
GO TO IBB 

IBS WRITE (POUT.B) , 

WRITE (POUT, 2) IDENT(t) 

IBB Site (POUT.ISI MENTIS). INPIIR>.IHI«>.TU.TCU.PL02.IR>. 
ITNETA(M).Pt(t.IR> 
ir (IR.EQ.NL) GO TO IBS 
IR«IR*I 
GO TO IBB 






non 


UMTC (POUT. ID IDCNTt*). (CODESt I . d. Ill) • d*l »2 ) >VrP( IR) . VVI>( IR 1 1 
itcoDESt2.d.iR).d.i.2).vsctiR).vrstiR).vvstiR) 
ir IIR.EO.NMS) CO TO izt 
00 TO tit 

READ tOt 01 MO CONDITIONS 

12t IR*t 

>} * »rIn*«PTM It) IDENT(Sl.NBStl.IR).NBSt2.IR).NBS().IA).THCS 
MRITCtPOUT.lt) IDENTt D.NBSt 1 .IR) * ,THCS 

■CRD (PlN.lt) IDCNTt D.HBSt l .IR) «HBSt Z.IR) .tttjt S.IR) 
limit (ROVT.lt) IDCNTt } ) .HBSt I . IR) .HBSt 2. IR) #t1BS( J. IR ) 
RCtDtPIH.lt) IDCNTt Il.tDBStl.IR). 1*1 ■*} 

MRITC tPOUT.lt) IDCNTt J).tDBS(I.IR). 1*1. t) 

IP (IA.EQ.HIC) 00 TO l*t 
00 TO lit 

It! CONTINUE 

RCRD OUTPUT OPTIONS 

RCRDtPIN.lt) IDCNT(S).IOUT 
MRITC (POUT. It) IDCNTt) ).IOUT 
Z FORMAT (IX.ltNIDCNTU) *.tt) 

i rotntT tzttt) 

• FORMAT (*' THERE IS R NI* UP IN THE LRYCR PROPERTIES CARD*) 
t FORMAT (At. Jit) 

It rotntT (At. Sit) 

12 rORNAT (Lt) 
u roRnAT tw) 

It rORNAT (At.2It.jrt.il 
u roRHAT (At.2it.SFti) 

It FORHAT (At.2AA«2rt.2.2At. JFR.Zl 
II rORNAT (At»2A*.2Ft.2.2A*.lFt .2) 

It rORNAT (At.TFl.t) 

It rORNAT (At.irt.A) 

2t FORMAT (At. It) 

21 FORMAT (At. It) 

22 FORMAT (IS) 

21 FORMAT ttCIt.l) 

RETURN 

END 


t 
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non no 


c „..5S!KS55l!5.«!S2!i5KiKiir;!;SKi^S5Si5?5i 

C SUMOUTINC TO SUPPLY VARIATIONS IM CONSTITUENT PROPERTIES 

g ....S5«555S555-5J-l25-S552-25-SS!:^-!:222!:i-5S2irSS2Si 

INTEGER PIN 
ORTA PIN/4?/ 

DINENSION°OUnH IBI.PrPttl .II.PPSIII . I I.PttPl U. I ) .PttS( U . 1 1 . 

Mli EACH^ PROPERTY HNICtt CORRESPONDS TO R BOOLEAN WITH VALUE 'TRUE* 

DO SB J • l.NTtS 

OENERATE TIBER PROPERTIES 

READ(PIN.IBBI) bool.shean.stdev 
IN. NOT. ROOD GO TO S 
CALL ORAHD(Xt) 

CALL URAND(X2) 

CALL NORMXI .M'SnCAN.STDEV.CrPM > 

PTPU.J) ■ EPPI1 

C S READIPIN.tBIII BOOL. SHEAR. STDEV 
IN. NOT. BOOL I 00 TO A 
CALL URANOIXll 

CALL mm(xi!xs.snERN,STDEV.crp>2) 

PTP(A.J) • EYP« 

C A READ! PIN. 1BBII D001.SHEAN.ST0EV 
in .NOT. BOOL) GO TO ? 

CALL GRAND! XI I 
CALL UAAHDIM) 

CALL N0RIMXI.X2.SnEAN.STDEV.GrPm 
PTPO.J) • 0PPI2 

C 7 READ! PIN. IBB II ROOL. SNEAN. STDEV 
IP(. NOT. ROOD GO TO B 
CALL URANOIXll 

call norhixi!x*.sheah,stdev. OTPII) 

PTP(t.JI • orm 

C R READ! PIN. I BB 1 1 ROOL. BETA. ALPHA 
IN .NOT. BOOL) GO TO f 
CALL URANDIRI) 

CALL WEIBIX). ALPHA. BETA. STPT) 

PFPI It. J) • STPT 

c 

I REAOI PIN. IBBI I BOOL. BETA. ALPHA 
IN .NOT. BOOL ) GO TO IB 
CALL URANDIRI) 

CALL HEIBIXI.ALPHA.DETA.SrPC) 
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prptis.J) • 

c 

It CONTINUE 

C OEMERATE MATRIX PROPERTIES 

C II READ (PIN. ltd ) BOOL • SHE AN ■ STDEV 
1F( .NOT. BOOL) 00 TO 21 
CALL URAHD(Xl) 

CALL NORM XI f XI . SNEAN . STOEW . EHNP ) 
pnp(j.j) « EtiHP 

C II BEAD(PXN.lMl) BOOL. BETA. ALPHA 
irl.NOT. BOOL) 00 TO 22 

CALL UBAND(Xl) . 

CALL UEIDIXl. ALPHA. BETA, SMTP) 
PflP(T.J) ■ SNTP 

C 22 READ(PIN.IRR) ) BOOL. BETA. ALPHA 
IP). NOT. BOOL) 00 TO IS 

CALL HEIBIXiIaLPNA.BETA.SHCP) 
PHP( IB.tl) ■ SNCP 

C IS BEAD(PIN.IBtl) BOOL. BETA. ALPHA 
irt.NOT. BOOL) 00 TO 14 
CALL URAND(Xl) 

CALL UEIB(X1. ALPHA. BETA.SHSP) 
PtlP(U.J) • SHSP 

C 24 CONTINUE 

C BEHIND PIN 
SB CONTINUE 

IBB1 r0RHAT( 14X.L4.2E2B.1B) 

RETURN 

END 



153 



APPENDIX B 



155 


This appendix outlines the theories and equations in the ICflN 
program that are used in this project. In the first section on 
conposite micronechanics, the elastic and thermal properties of a 
composite ply are defined with respect to its principal rater ial axes. 
The rext section, devoted to laminate theory, contains the 
transformations and summations of ply properties used to arrive at 
laminate properties. The last section contains a brief discussion of 
tte failure criteria. 

1. Conposite micronechanics 

Tte theory for calculation of the properties of a unidirectional 
fiber composite ply based on the properties, volume fractions, and 
orientation of its constituents is knovm as composite micronechanics. 
In this section, the subscripts f, m, v, and / represent fiber, matrix 
void, and laminate, respectively. The symbolic notation and the 
equations used are summarized below: 

Volume fractions: 

k + k + k = 1 

£ m v 

Longitudinal fbdulus: 

E m * "An * “» E » 


Transverse Modulus: 
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E I22 ~ E /33 


1 (1 - E |n /E £ 22^ 


Shear Mbduli: 


712 


1 - ^ - G„/G fl2 ) 


723 


1 “ ^ t 1 ” G ij» /G T23^ 


Poisson’s Ratios: 


“/12 * »J13 * “m * - *»> 


“/23 ” k f u f23 * k n 


712 


2v - 
m 


E m E ' 221 


Coefficients of theriral expansion 

a ., . + k [ (a E /E_, . ) 
/II n L 1 m m *11' 


- “ml 


•in 


1 + k (E /E_, . - 1) 
m k ra /II ' 


a /22 • “»(* ' *r > 


1 + VA11 

Vn * k .( E „ - E /n> 


a 33 “ a /22 


a /22 k / 
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2. Laminate Theory 

This section describes the methods vtiich are used to calculate the 
elastic properties of laminates from the properties, orientation, and 
distribution of individual laminae. The elastic properties are then 
used to predict the response of the laminate to external loads. The 
methods used to predict stresses in the laminae under application of 
external loads are also described. Failure loads can be predicted by 
using these methods; as described in a following section. 

a. Gereralized Hooke’s Law 

Tte stresses acting at a point in a solid can be represented by the 
stresses acting on the planes nornel to the coordinate directions, or 
equivalently, on the surfaces of an infmitesi.nl cube as shea, in Fig. 
B-l. The stresses (o, ) on each face are resolved into three 
components: one normal stress and t«, shearing stresses. The first 
subscript refers to the direction normal to the plane in tfiich the 
stress acts and the second subscript to the direction in %*uch the 
stress acts. The stress components show* on the faces of the cube are 
taken as positive and can be taken as the forces (per unit area) exerted 
by the material outside the cube upon the material inside. A stress 
component is positive if it acts in the positive direction on a positive 
face of the cube. Thus normal tensile stresses are positive, and normal 
compressive stresses are negative. Mine stress components must be used 

to define the state of stress at a point, namely <* u , O 33 , * 23 » 

„ - and a . There are nine corresponding strain 

*31* a 12’ *32’ *13’ *21* 
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co m ponent s , following the sane subscript convention. 

For bodies in tfiich each strain component is a linear function of 

all six stress conponents, the generalized Hooke’s Law can be expressed 

a. . = E . . f. , 

lj ljkl kl 

wtere E is a fourth order tensor of elastic constants. For nine 
i jkl 

stress conponents are! nine strain conponents, there nust be 81 elastic 

constants defining E i j kl - Certain reductions in the number of 

independent constants for an anisotropic body are due to synuetry 

properties of the tensor By considering nonent equilibrium about 

the center of the cube, it can be shovm that at any point ®23 = ^32* 

<? 31 = or 13 , and o 12 = a 21 - Thus, E. jkl is synnetric with respect to the 

first two indices. Second, because the strains are synnetric (that is, 

e = e ). E. .. . nust be synnetric with respect to the second tvo 
ij ji 7 ijkl 

indices. This reduces tte nunber of elastic constants to 36. Further 
reduction to the final 21 elastic constants for a general anisotropic 
neterial is accomplished by assuming the existence of a strain energy 


density function, such that 


with the property 


u = Ut-ij) 


au 

ds . . 
ij 


From the generalized Hooke’s Law, 


= E ijkl c kl 


Partial differentiation with respect to yields 
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_a_ [_3U_1 _ E 

*R! Kjl " ijkl 

Sirce tte order of partial differentiation is iimaterial, 


a 

f dU 1 

a 


df ki 

[as. 1 
1 ij J 

= de. . 

ij 

KJ 


and the subscripts can be interchanged to yield 



so that 


E ijkl ' ^lij 

Thus tte first pair of subscripts in can be interchanged with the 

secorrf pair without any change in the values. The number of elastic 


constants is thus rediced to 21. 


b. Lamina Constitutive Relation 

Several simplifications to the generalized Hooke’s Law can be made 
for tte special case of a thin orthotropic material, viiich approximates 
a unidirectional fiber conposite lamina. By considering the invariance 
of elastic properties under coordinate transformation for planes of 
symretry, the tensor E. jfcl can be reduced to the following nine 
constants: 

E U11 E 1122 E 1133 
E 1122 E 2222 *”2233 

E = E 1133 E 2233 E 3333 

i jkl 

0 

0 

0 

It is now convenient to make the following notation changes: 



16 ® 




vhere the coordinate axes coincide with the syimetry axis of the 
naterial. For laminae that are assuned sufficiently thin, the through 
the thickness stresses are zero. Thus — Oj = ®, for plane 

stress. It is apparent that ® 

The stress strain relations for a thin unidirectional lamina are 


written 
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a l 9 U 9 12 ® *1 

a 2 " 9 12 9 22 ® *2 

t 12 . ® ® ^66 • T y i2 

using tte tensorial strain y Y 12 instead o£ the engineering strain Y^. 
The Q terms are knovn as reduced stiffnesses, i.e. 

E» 

9 11 = C 11 = 1 “ » 12»21 

v l2 E 2 WaiEj 

9 12 = C 12 = 1 ~ U l2 W21 1 _ WlS v 21 

E 2 

9 22 = C 22 = 1 - Ui 2 U2i 

9 66 = T ^ C H " C 12^ = °12 

vftere E lf E*. «!*. v 21 , and G l2 are the ply elastic constants, neasured 
with respect to the natural material system. It may be noted that only 

four of these constants are independent. 

The stress- strain relation above shows that there is no coupling 
between tensile are! shear strains, as long as the applied stresses are 
coincident with the principal material directions. Hovever, coupling 
appears vhen a lamina is tested at arbitrary angles with respect to the 
principal material directions. The general form of the stress-strain 
relation for any angular orientation of a lamina is considered next, 
c. Stiffmss matrix transformations 

A lamina is loaded along a coordinate system x-y oriented at some 
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angle ® with respect to the principal material directions as shown in 
Pig, b- 2 . Since stress and strain are second order tensors, they are 

transformed by 

°1 

o 2 = [T] 

T 12 

and 

e i 

*2 
1 

l T V 12 

where [T] is the transformation natrix for plane stress and plane strain 

transformed by clockwise rotation about the (3,z) axes, gi'^n by 

cos 2 ® sin 2 ® 2 sin® cos® 

sin 2 ® cos 2 ® -2 sin® cos® 

-sin® cos® sin® cos® cos 2 ® - sin 2 ® 
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expressions for tte elastic properties referred to the x y coordinate 
system. 

d. Elastic properties of laminates 

A num ber of assumptions are trade in laminate theory to obtain 
theoretical predictions. These are: 

1 m lamina are perfectly bonded and do not slip relatiw to 

each other 

2. tie bond between the laminae is infinitesimally thin 

3. the laminate has the properties of a thin sheet 

Ttese assumptions allow the laminate to be treated as a thin 
elastic plate. The classical hypothesis of Kirchhoff is applied to 
derive tte strain distribution throixjhout the plate under external 
forces. Because tte laminate is conposed of laminae oriented in 
different directions with respect to each other, the stress-strain 
equation for each layer (k) is defined as 

^11 ^12 ^16 
*12 ^22 ^26 
^16 ?2S ^66 k 

Thus for a given strain distribution, the stress in each layer can be 
determired . The strain at any point in a laminate undergoing 
deformation mist be related to the displacements and curvatures of its 
midplare*. The discussion which follows assumes that the laminate is 
thin. Kirchhoff plate theory is used in this formulation. 

The deformation of an arbitrary section of a laminate is shown m 
Fig. B-3. It is assumed that lines straight and perpendicular to the 
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midplane before deforiration remain so after deformation. This is 
equivalent to neglecting transverse shearing deformations. Comparing 
Fig. B-4(b) with Fig. B-4(a), in Wiich the normals to the midplane 
retrain perpendicular after defcrrsetiCR, it is seen that the upper and 
lower surfaces of the plate must not shift their relatixra positions. It 
is obvious that the resistance of a thin plate to such defortration is 
large, much larger than its resistance to defoliations perpendicular to 
the midp lane. 

It is assuned that the point B at the midplane undergoes 
displacements u®, v®, and w® along the x, y, and z axes, respectitraly. 
The displacement u in the x direction of a point C located on the normal 
ABCD at a distance z from the mi dp lane is given by 


u = u® - za 


where a is the slope of the midplane in the x direction, 

dw® 

a = dx 

The last two equations can be used to obtain the displacement u of an 

arbitrary point at a distance z from the midplane as 

dw® 

u = u * ‘ 2 ~dT 


Similarly, 


v = v 0 


z 


dw® 

dy" 


Since the strains nor mal to the midplane are neglected (plane 


strain), tl* displacement w at any point is taken equal to the 
displacement w® at the midplane. The strains in terms of displacement u 


and v are 
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du duo 5 2 Wo 

£ x ~ dx = dx Z Ox 2 
dv dv 0 d 2 w 0 

€ y ~ dy ~ dy Z ay 5 

du flv du 0 dv 0 a 2 Wg 

xy dy dx dy dx dxdy 

In terns of midplans strains and plate curvatures, the strains in a 

laminate vary linearly through the thickness. 


s 


*° i 


k 

X 


X 


X 

s 

__ 

5° 

+ z 

k 

y 


y 


y 

y 


V° 


k 

xy J 


xy J 

V xy. 


where midplane strains are given by 


cO 


du 0 

c 

X 


dx - 

cO 


dv 0 

£ 

y 


ay" 

y° 


dUg dVg 

y xy J 


dy + dx 


ani the plate curvatures by 


*yj 



The stresses in any (k) lamina can be obtained by subst itut ing the 


previous equation into the stress strain equation 


a 


X 


a 

— 

y 


T 

xy 

k 


^11 ^12 ^16 
^12 ^22 ^26 
^16 ^26 ^66 




r° 1 


k 




X 


x 




F° 

+ z 

k 




y 


y 




y° 


k 


k 


xy J 


xy J 




166 


e. Laminate Stiffness ffetrix 

Classical laminate theory provides a method of calculating the 
resultant forces and moments per unit length acting on the laminate by 
integrate the stresses acting in each lamina through the thickness (h) 
of the laminate. Resultant forces are obtained by 


N = 
x 


f ' 

N = 

7 J - 

V = J- 


h/2 

-h/2 

h/2 


a dz 
x 


a dz 


h/2 

h/2 


t dz 
xy 


Tte ironent resultants are obtained by integration through the thickness 
of the corresponding moments of stresses about the midplane: 


vj 

r 1 

xy J 


h/2 

0* 3 

-h/2 

X 

h/2 

0 3 

-h/2 

y 

h/2 

T 

-h/2 

xy 


The units of N , N , N are force per unit length and tl*, M y , are 

x y 

moment per unit length. The sign conventions are show in Fig. B-5. 

Using tte resultant force and moment relations, a system is defined 
that is statically equivalent to the laminate stress system, but applied 



167 


at the midp lane. Thus, the external loading has been reduced to a 
system that does not contain the laminate thickness or * coordinate 


explicitly. 

For a laminate consisting of n laminae (Fig. B-6), the resultant 
force-ronent system acting at the midplane can be obtained by adding 


integrals representing the contribution of each layer by 


N 

K r * 

N 

y J -i 


dz = 


Si J Vj 



-h/2 


„ - 5 r* |«; U 

flJVi T 

l xy J k 


Using the expressions for the stresses in the k-th lamina derived 
earlier, and noting that the midplane strains and plate curvatures are 
constant not only within the lamina, but for all laminae, it is apparent 
that they can be taken outside the integral sign. The stiffness matrix 
[Q] is constant within a lamina so it also can be taken outside the 


integration to giw 
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M 

X 


n 

*11 

*12 

*16 

W 

y 

= 

l 

^12 

q 22 

*26 

w 

xy J 


k=l 

^16 

^26 

*66 




*• 1 

r \ . 

I z dz 


JK 

J V 1 


y 

y° 

k J 


xy J 


n 

*11 *12 *16 


P 1* 


k 

X 

r 

_ , __ 


\ 2. 



2 

*12 *22 *26 


z dz 
K 


k 

y 

k=l 

*16 *26 *66 

k 

Vi 


k 

xy J 


Three new matrices, A. B. ., and D. . , are defined, viiere 

X J X J X J 


n 

fl ij - } ( V* "v - V.) 

k=l 

B ij * ~r 2 ( V* (h £ ' *4-»> 

k=l 
1 ” 

B u ■ — 2 ( ^j )k (h * " h *- i) 


k=l 

These new matrices, A, B, and D, simplify the resultant force and moment 
relations, and are knovn as the extensional, coupling, and bending 
stiffness matrices, respectively. The total plate const itut iw equation 
is then 


H 


A 

B 



M 

= 

B 

D 


k 


It may be recalled that in an orthotropic lamina with arbitrary 
orientation the shear stress is coupled with the normal strain and the 


normal stresses are coupled with the shear strain. In general, a 
resultant shearing force on a laminated plate produces midplane normal 
strains in addition to the expected shearing strain. Similarly, a 
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resultant nornl force will induce shear strains in addition to midplane 
nor net 1 strains • 

The nonzero coupling matrix B in the plate constitutive equation 
explains the coupling between bending and extension of the laminated 
plate. Thus, normal and shear forces at the midplane induce not only 
midp lane deformations, (and hence, midplane strains) but also twisting 
arxi bending, producing plate curvatures. Similarly, resultant bending 
and twisting noments induce midplane strains. 

f. Lamina stresses and strains 

Tim aim of the analysis of a laminated conposite is to determine 
the stresses and strains in each of the laminae forming the laminate. 
These stresses and strains are used with failure criteria to predict the 
loads for failure initiation for a laminate. The failure criteria are 
discussed in the section devoted specifically to that purpose. 

The strains in a lamina caused by external loading are a function 
of laminate midplane strains and plate curvatures, as previously 
discussed. Once the lamina strains are known, lamina stresses can be 
found using the lamina stress-strain law. Thus, the starting point for 
calculating lamina stresses is the determination of laminate midplane 
strains arxi plate curvatures in terms of the applied loading. The plate 
constitutive equation given previously can be inverted to give the 
midp lane strains and plate curvatures explicitly in terms of the 
resultant external forces and ironents. The result of the inversion 

process is 
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A’ B’ 1 f N 1 | A f B* 1 | N | 

= c* D’ I I n | S I B' D' I [ H J 
where A’ , B’ , and D’ are simplified forms of the inwrsion process 

results, a«l are functions of the A, B, and D matrices of the original 
form of tte plate constitutive equation. 

It is row apparent that with these equations, an analysis of a 
laminate subjected to external forces and nonents can be conducted : 

1. calculate midplane strains and plate curvatures 




2. calculate lamina stresses in global (x-y) system 



3. calculate lamina stresses in natural (longitudinal and 
transverse to fiber) system. 



The strain variations in a lamina are calculated in an analagous 
manner. Tte stress-strain variation is compared with the allowable 
stresses anrf strains in each lamina. Thus the load at viiich failure is 
initiated in one of the lamina can be calculated, as long as a strength 
criteria exists in terms of the lamina natural axis system. The 
f ornulat ion of lamina failure criteria is discussed in the next section. 
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3. Strength Theories 

It is assumed that the strength of a laminate imist be related to 
the strengths of the individual laminae. A simple failure criteria 
consists of evaluating the lamina strengths in their principal material 
directions subject to induced stresses or strains at the boundaries of 
the lamina. In this context, it is assumed that the lamina and its 
constituents behave in a linear elastic manner to failure. The strength 
analysis described here assumes that the behavior of each lamina in an 
arbitrary laminate is the same as the behavior observed in the natural 
axis system when the lamina is part of any other laminate under the same 
stresses or strains. In other words, it is assumed that the strength 
criteria for a lamina in plane stress is valid for any orientation of 
the lamina in a laminate. In the ICAN program, the lamina strengths are 
calculated using the expressions given below. 

Longitudinal tension 

s /ht = s rr * k f + k m E m /E fii ) 

LongituJinal conpression: 

The longitudinal compressive strength must be computed on the basis 
of three different criteria: 

a. rule of mixtures 

S /11C = S fC <k f + 

b. delamination 

S /11C = ( 13 S /12 + S »C^ 
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c. fiber microbuck ling 
S /11C = " 


F 2°. 


i - k £ (l - G m /G / * 12 J 


Transverse tension 


S ,22T * S. t (F<CT/EOOH) 


Transverse conpression 

S 

Transverse shear 


- s / PEHOW 
/ 22C mC 


C( F 1 " 1 + V G fl2> F 2 G /12 S mBJ 
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G F, 

I7» 1 


FACT 


where F^ and F^ are given by 



H* variable EEXOK is introduced for convenience: 

denom = [i - v£J(i - E m /E f22 )] ^ 1 ♦¥>(?- 1) ♦ l / 3 (V> - l) 2 


there ? is given by 


F 


1 


E /22^ 1 ~ “ E B» /E f22^ 



P = 


1 
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Tte variable FACT is used to correlate the strengths of HIE and Kevlar 
fiber conposites with the experimentally observed valtES. Since neither 
of these fibers is used in this work, FACT takes the valve unity. 
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Fig. B.3- Bending geometry in the x-z plane. 



o. Ooftoct od bar without thoar b. OtlbcM bar with thoor do- 
doformotioo* formation* 

B. 4-Shearing force deformations on straight cross section. 
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